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Abstract: For gauge theories which admit a dual string description, we analyze deep inelastic 
scattering at strong 't Hooft coupling and high energy, in the vicinity of the unitarity limit. 
We discuss the onset of unitarity corrections and determine the saturation line which separates 
weak scattering from strong scattering in the parameter space of rapidity and photon virtuality. 
We discover that the approach towards unitarity proceeds through two different mechanisms, 
depending upon the photon virtuality Q^: single Pomeron exchange at relatively low and, 
respectively, multiple graviton exchanges at higher Q^. This implies that the total cross-section 
at high energy and large is dominated by diffractive processes. This is furthermore suggestive 
of a partonic description where all the partons have transverse momenta below the saturation 
momentum and occupation numbers of order one. 
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1. Introduction 

High energy, small angle, scattering has always been the subject of much attention, since the 
early days of elementary particle physics. While the phenomenological Regge theory, centered 
around the 'soft Pomeron', has met with considerable success in reproducing the experimental 
data for hadronic collisions, a fundamental understanding of this theory rooted in QCD is still 
lacking and seems to even transcend our present capabilities — as neither perturbative QCD, 
nor lattice gauge theory, can be applied in the relevant 'Regge kinematics' (high energy, but soft 
momentum transfer). 

In perturbative QCD, the resummation of the logarithmically enhanced diagrams leads to 
the BFKL, or 'hard', Pomeron [|l| which predicts a power-like growth with s for the gluon oc- 
cupation number and the scattering amplitudes. The growth is however faster than observed 
in experiments, and is also problematic from a conceptual point of view as it violates unitarity. 
Within the past decade there has been considerable activity and some important progress to- 
wards understanding the unitarization of the BFKL Pomeron within pQCD. These efforts have 
been mainly oriented in two directions: (i) the calculation [|| and proper implementation |^] of 
the next-to-leading order corrections to the linear BFKL equation, and (ii) the inclusion of the 
non-linear effects, like multiple scattering and gluon saturation, responsible for the unitarization 
of the scattering amplitudes @, ||, 0, . 

Whereas weak coupling techniques are undoubtedly useful for the high-energy problems 
which involve, at least, one hard momentum scale — so like deep inelastic scattering at small x 
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and large photon virtuality (s ^ S> Aqqj-,), or hadronic scattering with large momentum 
transfer (s ^ |t| ^ ^qcd) — ' other hand it remains unclear whether such approaches 

can shed any light on the 'soft' ^ ^qcd) aspects of Regge phenomenology. For instance, is 
there any relation between the soft and the hard Pomerons, and if so, then how to describe the 
transition between the two ? Such questions are particularly difficult because of our inability to 
explore the transitive region where the QCD coupling constant becomes strong. It is therefore 
interesting to study model field theories which are analytically accessible at strong coupling and 
which share some of the basic features of QCD at high energy — like the existence of a partonic 
description at weak coupling and a Regge-like behavior at high energy, for any value of the 
coupling. 

= 4 supersymmetric Yang-Mills (SYM) theory with 'color' gauge group SU(A'^c) and the 
't Hooft coupling constant A = g^N^ treated as a free parameter offers a particularly interesting 
theoretical laboratory for this purpose. At weak coupling, the high energy behavior is just like in 
perturbative QCD — in particular, to lowest order, the BFKL equation is formally identical in 
QCD and 7\A = 4 SYM ^ — while at strong couphng and in the large hmit the AdS/CFT 
correspondence maps the theory onto a weakly coupled superstring theory which naturally 
exhibits Regge behavior (see, e.g., |12, ^ 0, |l^, 16|). Such a change of behavior with increasing 
the coupling strength bears some resemblance to the transition from the hard Pomeron to the 
soft Pomeron in QCD, although some important differences persists. (For instance, when moving 
from weak to strong coupling, the Pomeron intercept decreases in QCD, while it increases in 
the = 4 SYM theory.) Thus one can hope that some useful insights for the corresponding 
QCD problems may be gleaned from a study of the Regge regime in = 4 SYM theory and its 
variants. 

In this paper we study deep inelastic scattering (DIS) in strongly coupled = 4 theory 
at small-x where x ~ Q^/s is the usual Bjorken variable. The basic formulation was given 
by Polchinski and Strassler |15| in the strict large- A'^c limit {Nc ^ oo at fixed energy), where 
however no unitarity issue arises: being suppressed by a factor l/N"^^ the elementary scattering 
amplitude is always smaller than one, and is dominated by the exchange of a single 'Pomeron'. 
In this strong-coupling context, the 'Pomeron' is the t-channel object exchanged at tree-level 
and at high energy in the dual string theory — a reggeized graviton propagating in the 
curved space-time that the string theory lives in and which is asymptotically AdS^ x . 

Thus, in order to study unitarization, one should rather consider the high-energy limit 
s — > oo at large, but fixed, values of A'c, and this is what we shall do in this paper. This is again 
similar to QCD (at large Nc)^ where the unitarity corrections are suppressed by inverse powers 
of Nc, hence a flexible definition of the 'large- A'^c limit' is necessary to observe such phenomena. 
The problem of unitarization in the context of string theory is notoriously difficult, even at weak 
(string) coupling g1 oc l/N^, because of the need to resum multi-loop string amplitudes to all 
orders, and here we shall not attempt to provide a full analysis in that sense (see however the 
approaches in Refs. |l^ Rather, our goal is more modest and also more pragmatic: 

Following again the example of weakly-coupled QCD — where it has been shown m, H, H, H] 
that one can follow the onset of unitarity corrections via an appropriate extrapolation of the 
results at weak scattering (in that case, the single BFKL-Pomeron exchange) — , we shall 
compute, at strong coupling, the saturation line which separates the weak-scattering from the 



- 2 - 



strong-scattering regions in the kinematical plane {x,Q^) for DIS. 

Specifically, while staying in the weak-scattering regime — which at strong coupling, like 
at weak coupling, corresponds to sufficiently large values of for a given value x <^ 1 — , 
we shall first identify those contributions to the scattering amplitude which become dominant 
when extrapolating towards the unitarity limit, and then deduce the saturation line as the curve 
Q1{x) along which the amplitude is constant and of 0{l). What at a first sight may look as 
a rather straightforward exercise in unitarizing the single-Pomeron exchange of Refs. []T5| , |T^ ], 
turns out to be quite subtle, because of several issues: 

First, the formulation of the unitarity constraint is rather non-trivial in the context of DIS, 
where the scattering is initiated by an off-shell, electromagnetic, current (actually, an 7?.-current 
in = 4 SYM). Within the context of pQCD, this issue is most conveniently solved by using 
the dipole factorization of DIS at small x |^, |2^, where one can unambiguously identify the 
on-shell hadronic state which scatters off the proton, and for which the unitarity bound can 
be properly formulated: the quark-antiquark excitation of the virtual photon, or 'color dipole' 
(see also Sect. 2 below, and Refs. [^, 27| for more detailed discussions). It turns out that a 
formally similar factorization holds at strong coupling as well, with the color dipole replaced 
by a vector-field fluctuation of the background metric (the '7^-boson'), which propagates on- 
shell before eventually scattering off the 'dilaton' target ||l^. This makes it natural to enforce 
the unitarity constraint at the level of the elementary scattering between the 7^-boson and the 
dilaton. This constraint will determine our saturation line. 

Second, unlike what happens in pQCD, where the single (BFKL) Pomeron exchange controls 
the dipole amplitude in the vicinity of the unitarity limit for all the relevant values of x and 
Q^, the situation at strong coupling appears to be more complex, and also more interesting: 
there, the single-Pomeron exchange represents the dominant contribution (when approaching 
the saturation line, once again) only for sufficiently low values of Q^, up to a critical value 
Ql oc exp{41niV2/^|^ Q'^>QI, a new mechanism takes up the leading role, namely 

the dijfractive scattering via the exchange of two, or more, elementary (massless) gravitons. 
Note that — > oo in the strict large- A'^c limit, which explains why this new mechanism has 
not been included in the analysis in Ref. |15]. Interestingly, the picture of DIS at high energy 
and > Ql that we shall arrive is quite close to that of Regge scattering in superstring theory 
in fiat space, as studied by Amati, Ciafaloni and Veneziano [12|. 

In fact, for a strongly-coupled gauge theory having a dual string description, the graviton 
exchange is naturally more effective, at high energy and large Q^, than the single-Pomeron one: 
the amplitude for one-graviton exchange grows very fast with the energy, as s ~ because 
the graviton has spin j = 2, and, moreover, its long-range non-locality can efficiently match 
the large separation in scales between the highly virtual 7?.-current and the dilaton. (Of course, 
being real, the single graviton exchange cannot contribute to the DIS structure functions, but 
multiple exchanges can do so, via diffractive final states where the outgoing fields are separated 
in rapidity.) The reason why the single-Pomeron exchange is nevertheless found to dominate at 
relatively low is because, in the AdS^ geometry, the massless graviton is turned off at very 
short separations (in the radial direction of AdS^)^ by the curvature of space-time. Or, by the 
holographic principle of AdS / CFT (which becomes manifest in the factorization of DIS alluded 
to above), small separations in AdS^ correspond to low virtualities on the gauge theory side. 
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This transition, from a single-Pomeron to multiple graviton exchanges, with increasing 
has important consequences on the shape of the saturation line, correspondingly obtained as 
the juxtaposition of two curves which smoothly match with each other at = Q^, and also 
on the behaviour of the DIS cross-section in the vicinity of this line. We thus find that, at 
large Q^>Q^, the structure function F2{x,Q'^) is strongly peaked at the saturation momentum 
— namely, F2 oc N^Q^ when < Ql{x), but F2 oc 1/{N'^Q'^) when ^ q2(^£^ _^ which 
is suggestive of a partonic interpretation with all the partons living at transverse momenta 
k±_ ^ Qs{x) and having occupation numbers of 0(1). This is further supported by an analysis 
of the 'sum-rule' J dxF2(x,(5^) which expresses the energy-momentum conservation: at large 
Q^ this sum-rule appears to be saturated by values of x near the saturation line, i.e., the 
largest values of x at which one can still find the partons. Needless to say, at strong coupling 
such a partonic language must be used with care, since the standard relations between structure 
functions and parton distributions hold, strictly speaking, only in perturbation theory. It is 
nevertheless interesting that the picture of parton saturation suggested by our analysis appears 
as a natural continuation towards strong coupling of the corresponding results as weak coupling, 
as obtained in pQCD |2|, g |, |, 0, |, |, |7|. 



This paper is organized as follows: In Section 2, we give a self-contained review of high- 
energy DIS in QCD at weak coupling, in the framework of dipole factorization. The focus is 
on the definition and the construction of the saturation line, with the purpose of facilitating 
the subsequent discussion at strong coupling, and also the comparison between the weak and, 
respectively, strong-coupling scenarios. In Section 3, we describe our setup of DIS at strong 
coupling following [^] and demonstrate a limitation of the single-Pomeron exchange approxi- 
mation, which fails to saturate the energy-momentum sum rule at large . This points out 
towards the importance of the massless graviton exchanges, which are then discussed at length in 
Section 4. On this occasion, we use and extend the 'Pomeron' approach initiated in Refs. [14, ^] 
and more systematically developed in |16], with the purpose of demonstrating the reemergence 
of the massless graviton propagator at large Q^. We then describe multiple graviton exchanges 
at a heuristic level, using in particular the results in Refs. |12, 17, |l^ to deduce the onset of 
unitarity corrections. We summarize our results for the saturation line in a 'phase diagram'. 
Fig. ^ for gauge theories at high energy and strong 't Hooft coupling, to be compared with the 
corresponding diagram at weak coupling in Fig. ||. Section 5 is devoted to a discussion of our 
results, in relation with a possible partonic interpretation. 



2. Saturation momentum and geometric scaling in perturbative QCD 

As mentioned in the Introduction, deep inelastic lepton-hadron scattering (DIS) at relatively 
high virtuality ^ -'^qcd exchanged photon and in the high-energy limit s ^ 

represents a rather clean laboratory for a theoretical study of unitarization in the framework 
of perturbative QCD — in particular, for the calculation of the saturation momentum. More 
precisely, in what follows we shall work in the 'small-x' regime of DIS, characterized by 

C>2 1 
as <Cl, X = — <Cl, and In - > 1 . (2.1) 

s X 

A 'lowest-order calculation' in this regime requires an all-order resummation of the perturbative 
contributions of order (a^ ln(l/x))"', with n > 1. This resummation can be performed by 
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solving appropriate evolution equations, constructed within pQCD, which are generally non- 
linear, with the non-linear terms accounting for the 'unitarity corrections' — gluon saturation 
and multiple scattering |5|, |7|, ^ 30|. However, a linear approximation — the BFKL equation 



1^, — holds within an intermediate range of energies, in which the non-linear effects 
remain parametrically small. The boundary of the validity region for this approximation in the 
kinematical plane {x,Q'^) is known as the 'saturation line'; this is the curve = Q'^{x) along 
which the scattering amplitude is constant and of order one. Remarkably, it turns out that 
the position of this line can be determined via calculations based on the BFKL approximation 
alone, without a detailed knowledge of the non-linear dynamics responsible for unitarization 



Q, 21, In what follows, we shall briefly review this calculation to 'leading-logarithmic 



accuracy' (LLA), i.e., at the level of the leading-order BFKL equation. (See also Ref. [31, ^] 
for more accurate analyses, including NLO effects.) 

To LLA, the cross-sections in QCD at high energy can be computed within a factorization 
scheme known as 'A^j—factorization' (see p3| , p6[), which is consistent with the non-locality of the 
high-energy evolution in transverse coordinates. When applied to DIS, this is most suggestively 
written as dipole factorization : the virtual photon 7* fluctuates into a quark-antiquark pair in 
a color-singlet state, or 'color dipole', which then scatters off the gluon field inside the target. 
More precisely, the DIS structure function F2 is computed as 



with (T-y*„ the total 7*p cross-section : 



a^*p{Q\x) = / dz /dV V |^'a(2,r)padipoie(r,x), (2.3) 

J a=T,L 



where the subscript a refers to the polarization of the virtual photon (transverse or longitudinal) , 
and ^aiz, r) is the wavefunction describing the dissociation of the virtual photon into a qq pair 
with transverse size r = |r| and where the quark (antiquark) takes away a fraction z (respectively, 
1 — z) of the photon energy. For massless quarks, one finds 



|M/i(z,r)|2 = 4J^Qh\l-zfKi{Qr), (2.4) 

where = z{l — z)Q'^ , and K\ are modified Bessel functions, and T = (A''caem/27r^) "^jej, 
with Cf the electric charge of the quark with flavor /. The factor of Nc comes from the sum 
over the color degrees of freedom of the quarks. The integration over r in Eq. (2.3) is effectively 



restricted to r < l/Q by the Bessel functions, which are exponentially suppressed for Qr S> 1. 

Furthermore, Eq. ( p.3p involves the total cross-section cTdipoieix, r) for the scattering between 
a color dipole with transverse size r and the hadronic target. For the present purposes, it is 
convenient to take the target itself as a collection of dipoles, succinctly referred to as the onium. 
In its own rest frame, the onium is a collection of dipoles with sizes r' < Rq distributed according 
to a density function n(r') and which are homogeneously distributed in impact parameter space 
within a disk of radius Rq. Note that Rq plays the role of an 'infrared cutoff' for this dipole 
distribution, in the sense that this is the maximally allowed dipole size. To ensure applicability 
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of perturbation theory, we shall assume that l/i?Q S> ^qcd (but this assumption is often relaxed 
in phenomenological studies of DIS, in which case Rq is identified with the proton radius). Also, 
it is convenient to choose ^ l/-^0' which case the integral in Eq. (2.3) is dominated by 
relatively small dipole sizes r (for the projectile dipole), such that r <^ Rq. 

To the accuracy of interest, the dipole-onium cross-section can be evaluated as 

o-dipolc(?^, t) = / -^aDD{r,r',T)n{r') , (2.5) 



where r = ln(l/x) is the 'rapidity' and ai:,j){r,r' ,t) is the total cross-section for the scattering 
between two dipoles with transverse sizes r and r' separated by a rapidity interval r. The dipole 
density is chosen as 



n{r') = l^j^J e{Ro-r'), (2.6) 

where the power A plays the role of an 'anomalous dimension': indeed, the standard situa- 
tion in perturbative QCD (at least for sufficiently small r') is A = 0, corresponding to the 
bremsstrahlung of small dipoles. The associated distribution is not normalizable (it exhibits a 
logarithmic divergence at r' — > 0), but this poses no problem for the calculation of the cross- 
section ( p. 51 ) since cri^Bif, f' , t) vanishes sufficiently fast when r' — > (see below). A normalizable 
hadronic state would have A > 0, in which case n(r'), once properly normalized, would have the 
interpretation of the probability density to find a dipole of size r' in the target. In what follows, 
we shall consider the cases where A is positive or zero. 

We shall soon compute the cross-section ( |2.5D in the BFKL approximation and then use 
the result, in association with the unitarity constraint, to deduce an equation for the saturation 
line. In preparation of this, it is useful to open a parenthesis and explain our specific use of 
the unitarity constraint for the present purposes. (A similar use will be made in the subsequent 
sections in the context of DIS at strong coupling.) 

The unitarity bound is most easily formulated in terms of the scattering amplitude T(r, 6, r) 
at fixed impact parameter h. It then reads T(r,b,T) < 1, where the upper bound T = 1 (the 
'black disk limit') describes a situation where the scattering occurs with probability one. Note 
that we use conventions in which the S'-matrix is written as S" = 1 — T, with T a real quantity 
at high energy. (That is, our T corresponds to the imaginary part of the conventionally-defined 
'scattering amplitude', which becomes predominantly imaginary at high energy.) However, the 
previous formulae, Eqs. (|2.3| )-( ^^ , involve directly the cross-sections, which are obtained from 
the respective amplitudes after integrating over b: crdipoie{r,T) = 2j d^6r(r, 6, r), etc. The 
unitarity constraint then becomes more subtle: with increasing energy, a cross-section can rise 
indefinitely, even after the 'black disk' limit has been achieved at central impact parameters, 
because of the radial expansion of the black disk. This expansion is however quite slow, loga- 
rithmic in the energy (to comply with Froissart bound), which is indeed much slower than the 
corresponding power-like increase of the amplitude in the BFKL regime (see below). Hence, it 
is possible to neglect the Froissart expansion (or, at least, treat this in an adiabatic approxi- 
mation) in a first study of unitarization. Then, the unitarity bound reads crDD(r, r', r) < 27rr^ 
(with r> = max(r, r')) for the dipole-dipole scattering, and (Tdipoic{x,r) < 27ri?Q for the dipole- 
onium one. It is in fact more suggestive to introduce scattering amplitudes averaged over b, e.g., 
T{r,T) = crdipoie(r, r)/(27ri?Q), in terms of which the unitarity bound reads simply T < 1. 
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We can now explain our strategy for computing the saturation line: Starting in the regime 
where the scattering is weak, T ^ 1, we compute the (average) scattering amplitude from 
the solution to the BFKL equation and then extrapolate the result towards the unitarity limit 
T ~ 0{1). The saturation momentum Qs{t) is then obtained via the condition 

r(r,r) ~ 1 for r~l/Qs(r). (2.7) 

That is, 1/Qs{t) is the critical dipole size for the onset of unitarity corrections (in the form of 



multiple scattering) in dipole-target scattering at rapidity r |2C, 21 1. This is also the critical 
momentum scale for parton saturation in the target wavefunction [ p!9| , 28, 29, ^, |9|, in a frame 
where the target has rapidity r : for transverse momenta k± ^ Qs{t), the gluon occupation num- 
bers saturate at a value Ug ~ l/(asNc), which in turn implies that the (anti)quark occupation 
numbers saturate at a value Ug ~ 1. (Recall that, at small x and relatively low k±, the quark 
distribution is driven by gluons, which dominate the small-x part of the wavefunction.) This 
profound link between dipole unitarization and gluon saturation in pQCD can be understood as 
a consequence of the relation T(r, r) ~ asUg^r, k± ~ l/r), valid in the BFKL regime, between 
the dipole amplitude and the gluon occupation number. 



By using Eq. (2.7) together with the BFKL solutions to be shortly presented, we shall 
compute the saturation lines for both dipole-dipole, and dipole-onium, scattering. For more 
clarity, let us anticipate here the main conclusions that we shall arrive at: (i) The value of the 
saturation momentum depends upon the nature of the target, but its evolution with increasing 
energy does not — the evolution is universal, (ii) The saturation momentum of the onium 
coincides with that of its 'softest' component — the constituent dipole with maximal size r' ^ Rq. 

We now close the parenthesis and return to the calculation of the cross-section in the 
BFKL approximation. When applied to the dipole-dipole cross-section aDD{r,r' ,t)^ the BFKL 
equation [|| reads (with ag = asNc/ir) 



dr 



o-DD(r,r ,r) = — d z 



{ - a-DY){r,r' ,t) + 0-00(2, r',r) + (TDD(k - z\,r' ,t)] . (2.^ 



This equation has been written here as an evolution in the projectile size r for a given target 
size r' . (Of course, the function aDD{r,r' ,t) is symmetric in r and r'.) It has a simple physical 
interpretation: in one evolution step (r ^ r + dr), the dipole r can split into two dipoles, 
with sizes z and respectively r — z, which then interact with the target. The 'BFKL kernel' 
^ [r'^/z'^[r — z)'^] represents the differential probability for such a splitting to occur per unit 



rapidity. To solve Eq. (2^), we also need an initial condition at low energy (r = 0). To the 



present accuracy, this is computed as the exchange of two gluons, which yields 

aoD(r, r', r = 0) = 27ray^ (^1 + In , (2.9) 

where r< = min(r, r') and r> = max(r, r'). This cross-section shows 'color transparency' — it 
vanishes like when r< ^ — , meaning that a very small dipole scatters only weakly. By 
using Eq. (^.g|), one can easily check that, for r = and r <^ Rq, the dipole-onium cross-section 
(|2.5|) is dominated by relatively large target dipoles, with r <^ r' < Rq. (When A = 0, the large 
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dipoles dominates over the small (r' < r) ones by a large logarithm lii{RQ/r).) As we shall later 
see, the feature is preserved by the evolution with increasing r. 

The solution to the BFKL equation with the initial condition ( |2.9D is most conveniently 
obtained in Mellin space, where it reads 

c 

where x{l) is the eigenvalue of the BFKL kernel in Mellin space, generally referred to as the 
BFKL characteristic function 

x(7)=a.{2V'(l)-^(7)-^(l-7)}, ^^(7) = dlnr(7)/d7, (2.11) 

and the integration contour C runs parallel to the imaginary axis with < Re(7) < 1. Note 
the symmetry property x{l) = x(l ~ 7); which reflects the conformal invariance of the BFKL 
evolution. In view of this, one can check that Eq. ( p.lOj) is symmetric under r r', as it should. 
The inverse Mellin transform in Eq. ( p. 10 ) can be evaluated via a saddle point approxima- 



tion, in which the value of the saddle point depends upon the balance between the kinematical 
variables a^r and p — p' = \Q.{r''^ / r"^) . Here, we have introduced the logarithmic variables 
p = ln(i2Q/r^) and p' = \\i{B^/r''^) which are convenient for what follows. At this point, the 
reference scale inside the logarithms is arbitrary (it anyway cancels out in the difference p — p'\ 
but for later convenience we have chosen it as the onium size i?o- 

Consider first the formal high-energy limit olsT 00 at fixed p — p' (this is formal since 
the respective solution violates the unitarity bound already for relatively low energies, as we 
shall shortly see). Then the saddle point 70 is determined by the minimum of the characteristic 
function : 

x'(7o) = =^ 70 = 1/2, (2.12) 

and Eq. (|1^) can be evaluated by expanding x(7) around 70, to second order ('diffusion 
approximation'): writing 

7 = l/2-ii/, (2.13) 

one obtains (notice that x(l/2 — if) is an even function of z^, by conformal symmetry) 

xiXll - iv) = ujQ - D^v^ ^ Oiy"-) with a;o = 41n2a^, L>o = 14C(3)a, , (2.14) 

and then the Gaussian integration over v can be easily performed to yield 

. , ^ 2 e^o^ , 5 ,5,1/2 f ln2(rVr'2)l 
(^DD(r, r',T) ~ Sag ^ (rV'2)'' exp<^ '- \ . (2.15) 



This expression exhibits the 'BFKL Pomeron', i.e., an exponential increase with r at fixed r 
and r'. Clearly, when extrapolated at large r, this behaviour violates the unitarity bound. 
For instance, when r ~ r', this happens at a critical rapidity Tcr ~ (I/ojq) ln(l/ag), which is 
parametrically not that high. (Recall that the 'high-energy' regime starts only at r ~ l/a^.) 
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It is nevertheless possible to explore higher rapidities r > r^r while staying within the limits 
of the BFKL approximation provided we evolve along an oblique direction in the kinematical 
plane p — t (see also Fig. |^) : when increasing r, we should simultaneously increase p (i.e., 
decrease the ratio r/r'), in such a way that the scattering amplitude remains small: T(r, r', r) <C 
1. In particular, the saturation line p = Ps{t) corresponds to a direction of evolution along which 
the amplitude is constant and of order one, cf. Eq. (^). 

This also means that, in order to follow the saturation line, one must abandon the symmetry 
between the projectile and the target dipoles (r r'), which has been manifest so far. If the 
target dipole size is fixed to r', then with increasing r we must evolve towards smaller and 
smaller values for r. It is then natural to write (Tdd('^5 ^^'^ ''") = 27rr'^r(r, r', r), where r < r' and 
the unitarity bound corresponds to T{r,r',T) = 1. The BFKL approximation holds so long as 
T{r,r',T) <^ 1, but it allows us to approach the saturation line from small r (large p) values, 
and thus compute the saturation momentum according to Eq. ( p.T] ). To that aim, however, 
one cannot rely on the Pomeron saddle point, Eq. ( p.l2| ), anymore, rather one must determine 
the saddle point 7^ corresponding to an evolution along the saturation line (p = Ps{t)). By 
combining the saddle point condition: 

Tx'{ls) = -(.Ps-p'), (2.16) 



with the condition that the amplitude (2.10) be approximately constant along the saturation 
line: 

Tx(7s) - (1 - ls){Ps - p) - H^/al) = 0, (2.17) 
(we have only kept the dominant parametric dependencies upon r, ps and a^), one obtains 
X'ils) _ 1 



Xils) 1 - Is 



7, ^ 0.372, 



( '\ '1 ln(l/a2) _ x{ls) _.CQQ- ^^18^ 
Ps{T,r) = p +VsT with Vs = 4.883qs . (2-18) 

1 - 7s 1 - 7s 

(We have assumed here that r, and hence ps, are large enough for the term ln(l/o!^) to be 
treated as a small perturbation.) Note that 7^ is a pure number, independent of either r or cxg, 
corresponding to the fact that /9s(t) is a straight line, with slope Vg ~ 0{6is)- 

For p larger than p^, but not much larger, one can estimate the amplitude by expanding 
Eq. ( 2.10| ) around 7^ : writing 7 = 7s — iv, and expanding to second order in v, one finds 



Txil) - (1 - 1)P - ln(l/a2) ~ -(1 - 7,)(p - p,) - iv{p - p,) - D,tu\ (2.19) 

where Dg = x"(7s)/2 ~ 24.26as- This limited expansion is valid so long as 1 < p — ps <^ DgT. 
In this range, the amplitude is obtained by performing the Gaussian integration over i/, as 

where the saturation momentum (here, for a target dipole with size r') 

Ql{r\T) = (a2)V(i-7.)_le-- (2.21) 
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appears as the natural reference scale for measuring the projectile dipole size. When viewed as a 
function of rQs, the amplitude near saturation, Eq. ( |2.20| ), is universal, i.e., independent of the 
specific properties of the target, but uniquely fixed by the evolution. Note that, in Eq. ( |2.20| ), 
7s plays the role of an anomalous dimension. The evolution of the saturation momentum with 



increasing r, cf. Eq. ( |2.2l| ), is universal as well, and controlled by the 'saturation exponent' 
Vg- The above calculation yields Vs ~ C(l) for as = 0.2 0.3, a rather large value that would 
be inconsistent with the phenomenology at HERA or RHIC. However, it turns out that, after 
taking into account the next-to-leading order BFKL corrections ||2[ |3|, this value is reduced to 



Vs ~ 0.3 |31], which is indeed in the ballpark of the various phenomenological analyses (see, e.g., 
H, P, H for studies at HERA). 



We now return to the problem of dipole-onium scattering, cf. Eqs. (2.5)-( |2l^ ), with the 
purpose of computing the corresponding saturation momentum Qs{t). As before, we shall 
approach the saturation line from the weak scattering regime, i.e., from relatively small dipole 
sizes r <^ 1/Qs{t). This procedure meets with a subtle point though: although the overall 
scattering is, by assumption, weak, T{r,T) <^ 1, this is not necessarily so for all the individual 
dipoles which compose the onium and which contribute to the convolution in Eq. ( |2.5| ). The 
problem comes from the relatively small dipoles which, according to Eq. ( 2.211) , develop large 



saturation momenta (meaning that they evolve into 'hot spots' with high gluon density), off 
which the projectile dipole can strongly scatter. Hence, when evaluating the contribution of 
such small dipoles to the total cross-section, one cannot rely on the BFKL approximation 
anymore. Still, as we shall shortly argue, this brings no serious complication, since, after being 
properly unitarized, the small dipoles give negligible contributions to the overall cross-section 
in the vicinity of the saturation line. 

To see this, one can use the following, piecewise, approximation for the dipole-dipole cross- 
section at high energy. 



(^DDir,r',T) « 27rr> < 



''T{r^Qsir>,T), t) for r< «; 1/Qs(r>,r) 

(2.22) 

1 for r< > 1/Qs(r>,r) . 



In the first line, T{r<^ Qs{ry,T),T) is the amplitude in the weak scattering regime, as given 
by Eqs. ( 2.2C| )- (|2.2l| ) after replacing r — > r< and r' — > r>. The second line represents the 



black disk regime. The precise transition between these two regimes is not under control. 



but this is irrelevant for the present purposes. By inserting Eq. ( 2.22 ) into Eq. (|2.5| ) and 
performing the integration over r', one can estimate the dipole-onium scattering amplitude 
r(r, r) = (Tdipoic(?'5 ''")/2vri?Q. The very small target dipoles with r' < r are easily seen to be 
irrelevant, since their contribution to T{r,T) is of order r^/B^ <ti 1. (Recall that we are only 
interested in contributions which can approach the unitarity limit T = 1 when r 1/Qs{t).) 
To evaluate the contribution of the larger dipoles with r < r' < i?0) it is again convenient to use 
logarithmic variables {p = ln(i?Q/r^), etc.), and notice that Eq. ( |2.18| ) can be rewritten as 

Ps{T,r') = p' + Pso{r), Psoir) ^ VsT - ^"/^^"'^ , (2.23) 

where PsoiT) represents the saturation momentum for a dipole with the maximal size r' = Rq. 
We anticipate here that in the regime of interest we have p > Pso{t). Then the contribution of 
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the target dipoles with p' < p (i.e., r' > r) can be evaluated as 

T{p,T)^ J dp'e~^'-''^')^P-p'~'>^^h- '''"^^T' ^-(1+^)P' + I dp'e-(i+^)''', (2.24) 

p-pso 

where the first (second) term in the r.h.s. corresponds to large dipoles which scatters only 
weakly (respectively, to relatively small dipoles for which the unitarity limit has been reached). 
The first integral in the r.h.s. involves the overall exponential factor 

e(i-7.)p'e-(i+^)p' = e-(^+^^)''' (2.25) 

where the exponent A + 7^ is strictly positive in pQCD (recall that A > and 7^ ~ 0.37). It 
is then easy to check that, so long as p — pso < 2(A + 7s)Z)sT, this integral is dominated by its 
lower limit p' = 0, i.e., by target dipoles of the largest possible size r' ~ Rq. As for the second 



integral in Eq. ( 2.24 ), this is dominated by its own lower limit, at p' = p — pso- One thus finds 



T(p, t) — e-(i-^=)(^-^^o) e-^^^ + ^-{i+A)ip-p^o) . (2.26) 

^ ^ A + 7, A + 1 ^ ^ 

As just mentioned, A + 7^ > 0, hence the second term in the r.h.s. is exponentially suppressed 
w.r.t. the first one. Thus, finally (up to a slowly varying prefactor), 

T{p, t) ~ e-(^-'fs)ip~Pso) e" 'iofr for p - Pso < 2(A + Js)DsT , (2.27) 

which is essentially the same result as for a target made with a single dipole of size Rq (compare 
to Eq. ( p. 20 )). In particular, the saturation momentum of the onium coincides with that of its 



largest dipole component (cf. Eq. ( 2.21| )) : /Os(t) = Psoi'T'), or 



Q'sir) = (a2)i/(i-7.)_Le--. (2.28) 

This is quite intuitive: the onium starts to look 'black' as a whole only when a large dipole, 
which covers all (or most) of the hadron disk, has evolved into a system with high gluon density 
on the resolution scale of the projectile. Smaller dipoles reach the unitarity limit much faster (on 
that particular resolution scale), but their contributions to the total cross-section are suppressed 
by their small area^. 

The high-energy dipole amplitude in Eqs. ( 2.20 ) or ( |2.27| ) exhibits an interesting struc- 



ture, which can be better appreciated by comparison with the low-energy amplitude (|2.9| ), or 
the (purely) BFKL prediction at high energy, Eq. ( 2.15| ): Besides the power-law r^^^""^') with 



anomalous dimension 7^ « 0.37, there is also a Gaussian factor describing diffusion in the log- 
arithmic variable p ln(l/r^). This shows that, when increasing r, the amplitude gets built 
via a diffusive process which gets most of its support from the region < p — < \/4:DsT 
above the saturation line. For r large enough, such that Q'^{t) S> ^qcb^ ^^^^ region lies fully 
in the perturbative domain, thus confirming the internal consistency of our present approach. 
The physics of saturation/unitarization eliminates the symmetric diffusion characteristic of the 
(pure) BFKL evolution, cf. Eq. ( p. 151 ), which would invalidate the use of perturbation theory. 

^Of course, the physical picture would be different if A was negative and large enough, in such a way to bias 



the dipole distribution (2.6) towards small sizes. Then, the onium could become black as a whole (on a given 



resolution scale) by getting covered with many small dipoles which are individually black. 
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Figure 1: A 'phase diagram' for the high energy evolution in pQCD in the presence of unitarity correc- 
tions and gluon saturation. (We recall that as = asNc/ir.) 



In particular, within the more restricted window p — <^ ^JADgT (close to the saturation 
line), the diffusion term in Eq. ( p. 27 ) can be ignored, and then the amplitude shows geometric 
scaling, i.e., it depends upon r and r only via the single variable rQs^r) : 



r(r,r) 



-(l-7s)(p-Ps) 



for 



T . 



(2.29) 



This scaling property is an important consequence of the non-linear dynamics responsible for 
unitarization on the shape of the amplitude in the weak scattering regime 'above saturation' 
{p > ps). Via the dipole factorization (^I^), this property gets transmitted to the DIS cross- 
section in the high regime at > Q1{t) : 



i^2(Q^r) 



for QliT)<Q^<^QliT), (2.30) 



where Q'^{t) = Qg(r)e^^-^''^ is the upper bound of the scahng window, cf. Eq. (|2.2g| ). 

It is furthermore interesting to note the behaviour of F2 for very small values of Q^, deeply at 
saturation (Q^ ^ Q'^{t)). Then, the virtual photon cross-section is dominated by large dipoles, 
with 1/Q^(r) < < l/Q'^, for which the black disk limit is saturated: o"(iipoio(^) ''") ~ 2ttRq. 
One then finds (the logarithm is generated by the integral over r^) 



"2 l^cKoln—^ 



for 



Q'«Q^(t), 



(2.31) 



which exhibits geometric scaling too. Remarkably, it turns out that a similar scaling has been 
observed indeed ISq] in the HERA data for DIS at small x < 0.01 and up to relatively large 
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values of ~ 100 GeV^ (well above the saturation momentum at HERA, estimated in the 
ballpark of 1 GeV). The 'phase-diagram' for high energy evolution in pQCD which emerges 
from the previous considerations is graphically summarized in Fig. |^. 

3. DIS at strong coupling 
3.1 The setup 

With this section we begin the study of the most interesting problem for us here, namely DIS 
at strong 't Hooft coupling and high energy. Specifically, our physical gauge theory will be a 
deformation of the = 4 SYM theory in the regime where the gauge coupling is weak, g <C 1, 
but the number of colors is large, Nc ^ 1, ui such a way that the 't Hooft coupling A = g^Nc 
(the relevant coupling for perturbation theory at large A^c) is large : A ^ 1. The 'deformation' 
refers to the introduction of an effective infrared cutoff which breaks down conformal symmetry 
and mimics confinement. 

According to the AdS/CFT correspondence [^], the A^ = 4 SYM theory is admitted to 
have a dual description in terms of the Type HB string theory living in the 10-dimensional 
AdS^ X space-time, with the following metric: 

ds^ = -^{-dxl + Axl + dxl + dxi) + -^du^ + R^d^l . (3.1) 

In this equation, xq and x = {xi,X2,xs) are the time and, respectively, three spatial dimensions 
of 'our' physical Minkowski world, u is the AdS^ radial coordinate (0 < tt < oo), to be below 
referred as 'the 5th dimension', and dO| is the angular measure on S^. Furthermore, R is the 
common radius of AdS^ and S^, and it is related to the string length Ig via the 't Hooft coupling 
of the dual gauge field theory: R = /^A* = a'^X^, with a' = l1 the Regge slope. The physical, 
4-dimensional, space is identified with the boundary of AdS^ at u ^ oo, where the AdS^-m.etv\c 
becomes conformal to the standard Minkowski metric rj^'^ . The string coupling constant gg and 
the gauge coupling g are related to each other as gs = g'^ / 4:1t = X/AttNc, which shows that in the 
'strong coupling' regime of interest here (A ^ 1 but g <^ 1), the string theory is weakly coupled 
and hence it can be treated in a semi-classical approximation. 

We have previously mentioned the 'deformation' of the AA = 4 SYM theory which is neces- 
sary to mimic confinement. In the dual string theory, this translates into a modification of the 
metric at small u ^ uq = AR^, where A is a parameter which sets the scale of light hadrons 
('glueballs'). The precise way of modifying the geometry in the infrared is an issue of debate 
which we do not enter. (As we shall see, the dominant behaviour at high energy is insensitive 
to such details.) In what follows, we shall simply restrict the 5th coordinate to values u > uq 
(in practice, we shall be mostly interested in n ^ uq), and assume that the spectrum of the 
supergravity theory includes states which are localized near u = uq, corresponding to massive 
'glueballs' on the gauge theory side. Note that, because of the logarithmic nature of the metric 
in u (ds^ oc du^/u^), massive 'supergravity' states exhibit power-law decays at large u, to be 
better characterized later on. 



To closely mimic the setup of DIS, we follow |15] and introduce a 'photon' as the gauge boson 



associated with a gauged U{1) subgroup of the 5'[/(4) 7^-symmetry. As for the target, we employ 
a 'supergravity' glueball of mass ~ A which originates from the Kaluza-Klein decomposition 
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of the 10 dimensional massless dilaton after compactification on S^. Still as in Ref. [^], we 
shall use indices M,N,... to denote all ten space-time dimensions, separating into fijU,... on 
the Minkowski boundary at ti — > oo (our physical space), m,n,... on AdS^, and a,b,... on 
Indices are raised or lowered in the standard way, using the metric Gmn which can be read 
off Eq. ( |3.l| ) : ds^ = GMNdx^dx^ . It is understood that for quantities living in the four- 
dimensional gauge theory (and hence also on the Minkowski boundary of AdS^), the tensor 
operations are performed with the Minkowski metric t?^*^ = (—1, 1, 1, 1). 

As standard in DIS, the structure functions are obtained from the imaginary part of the 
forward Compton scattering amplitude 



I J d'^xe'''-''{P\T{Jf'{x)r{0)}\P), (3.2) 



where J^(a;) is the density of the 7^-current and \P) denotes a (normalizable) glueball state 
with 4-momentum P^. The hadronic matrix element under the integral will be computed from 
semi-classical string theory, via the AdS/CFT correspondence. The recipe is as follows : 

The 7^-current at the Minkowski boundary (n oo) excites vector-like metric fluctuations 
in the 'bulk' (i.e., at finite u > uq) : 

u)Va{^), (3.3) 

where Q, denotes the angular coordinates on S^, Am{Xfj_,u) is a gauge boson field, and Va is 
the Killing vector on corresponding to the gauged U{1) subgroup. At n — > oo, the field 
^^(2:^,00) can be identified as the physical photon which couples to the 7^-current. This will 
be taken in the form A^(Xf^,oo) = n^e*'''^, with space-like four-vector : = q^q^ > 0. 
But inside the bulk, Am{x^,u) should be better viewed as a gravitational wave, as it describes 
a fiuctuation in the metric tensor. For definiteness, we shall refer to this field as the '7^-boson'. 
This 7^-boson propagates into the bulk until it scatters off the dilaton target — a supergravity 
state with wavefunction 

'f(x^,u,f]) = e*^-^^>(u,Jl). (3.4) 

Note that this is a plane wave with four-momentum P'^ in the four physical dimensions. The 
Bjorken-x variable for DIS is defined in the usual way: x = Q2/(2P . q). Still as usual, the 
high-energy limit (P — > 00 at fixed Q^) corresponds to small values for x {x <^ 1). 

As explained in Refs. |15], the Regge behavior starts to be seen when 1/x S> \/A where 



the energy is high enough to create excited string intermediate states. To leading order at large 
Nc, the string scattering can be computed at tree level, as the imaginary part of the superstring 
Virasoro-Shapiro amplitude for graviton-dilaton scattering (originally calculated in flat 
space, and heuristically extended to curved space in [|l5|) convoluted with the supergravity 
wavefunctions for the 7^-boson (Am) and the dilaton (<&). The result of can be compactly 
written as 

F,{x,Q') = ca'^^ I ^^jKi{QR'/u) + Kf{QR'/u)) 0) ~^^^Hu), 

(3.5) 
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where we have absorbed in c a constant coming from the integration over S^, as well as powers 
of A. (For simplicity, we do not keep track of powers of A in the prefactor. This is acceptable 
since we shall deal with an exponentially large factor e^. See below for details.) 

Let us discuss the various factors under the integral in Eq. ( |3.5| ) one after the other: 

(i) The modified Bessel functions i^o,! are a part of the '7?.-boson wavefunction'; that is, 
they arise via the solution Am{Xfj_, u) to the supergravity equation of motion for the metric per- 
turbation induced by the 7?.-current. This equation is Maxwell equation on AdS^: D"^Fmn = 0, 
with Fjnn = dm An — dn Am and D"^ the covariant derivative in the AdS^ background. {Am{x^,u) 
is the non-normalizable solution to this equation corresponding to the boundary condition 
Afj_{Xfj,, u ^ oo) = n^e*"^'^.) For u smaller than 

Uc ^ Qi^^ (3.6) 

the Bessel functions decays exponentially, meaning that a 'photon' with large cannot pene- 
trate deeply inside the bulk: the larger is, the closer the photon remains near the boundary. 
Since, on the other hand, the inverse power of u manifest in Eq. ( |3.5| ) favors small values of u, 
it is clear that the integral there is controlled by n ~ Uc. 

Remarkably, the 7^-boson wavefunction involves the same Bessel functions as the wavefunc- 
tion describing the qq excitation of a virtual photon in pQCD (cf. Eq. (|2.4| )). This suggests a 
kind of 'duality' between the dipole size r in pQCD and the 5th dimension u in AdS^, via the 
correspondence r *^ R'^/u. We shall identify further aspects of this correspondence later on. 

(ii) The dilaton wavefunction $(«) is peaked around u ~ uq and decays at large u like 
<I> ~ u~^, where A is the conformal dimension of the operator which creates this state. For 
glueballs corresponding to supergravity modes, A is a number of 0(1) which is strictly larger 
than 2. In the formal analogy to pQCD, the probability density |<I>(n)p ~ u~^^ corresponds to 
the density n(r) of dipoles inside the onium, cf. Eq. ( |2.6| ). Note that the overlap between the 
7^-boson and the dilaton wavefunction scales like u~'^^ ~ {l/Q"^)^, and hence it rapidly dies 
away when increasing Q^. 

(iii) The last comment shows that DIS would be strongly suppressed at large if there 
was not for the non-locality hidden in the middle factor, (l/x)^^^°'/^)^^, which describes the 
imaginary part associated with one string exchange in the curved, AdS^, space. In less formal 
notations, the action of this operator can be written as (with r = ln(l/x)) 

-] «>t$(n) = / dn'(de^|n')|$(n')p. (3.7) 

X J X J 

Notice the dominant energy behavior F2{x) ~ 1/x aX small x: this is the hallmark of the graviton 
exchange in the i-channel. (At high energy, the exchanged string is dominated by its graviton 
component; see Sect. 4.2.) Furthermore, the operator A2 in the exponent is the Laplacian 
describing diffusion in the 5th dimension, for the spin-2 exchanged graviton. This operator 
appears because of the finite momentum transfer in the u-direction, although the momentum 



transfer in the physical four dimensions is strictly zero [15, |l^. Explicitly, this operator reads 

(I) =7?'"' 



A.-(^rvSm" = 4(9^1) (3.8) 
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where Vq is the scalar Laplacian on AdS^ and we have defined p = ln{u^ /uq). By also using the 
integral representation of the 'heat kernel' exp{5p}, one arrives at the following representation 
for the expression in Eq. (| 



dp' e^T^ e-^''', (3.9) 



2VttDt Jo 

where we have used /a' = \fX and introduced 

ujo = l-^, and ^ = ^- (3.10) 

Eq. ( |3.9D result shows a Pomeron-like increase with r, with an intercept 1 + = 2 — 2/\/A 
which is slightly shifted from 2 due to the curvature of AdS^ [|14|, |l^, |l^ . 

One should notice the formal similarity between the stringy 'Pomeron' amplitude in Eq. ( |3.9| ) 
and the cross-section ( p.l5| ) which describes the BFKL Pomeron in pQCD. In this analogy, the 
fifth coordinate u plays, once again, the role of the inverse dipole size 1/r. Like for the BFKL 
Pomeron, the non-locality associated with the single Pomeron exchange in AdS^ is diffusive in 
p ~ Inn-^, with the rapidity r playing the role of the evolution 'time'; accordingly, its effect 
becomes important only for very high energies, exponentially large in the relevant coupling. 
Namely, the diffusive radius Dt in Eq. ( |3.9D becomes of 0(1) when r ~ \/A, or 1/x ~ e^ . For 
comparison, the corresponding estimate in pQCD reads r ~ l/cis, or 1/x ~ e^/"*, which at a 
first sight may look as an astronomically high energy, but in fact is is not: a<j is not that small 
in real world {oig = 0.2 0.4 for the high-energy experiments), and rapidities of order l/oig or 
even larger are within the reach of the present-day accelerators. 

Up to an overall factor of 0{1/N^) that we shall shortly comment on, Eq. ( |3.9| ) can be 
viewed as the string theory analog of the dipole-onium scattering amplitude in pQCD, in the 
BFKL approximation. The dipole (a qq excitation of the virtual photon with a given transverse 
size r) is now replaced by the 7?.-boson Am{u) (the metric fluctuation induced by the 7^-current 
at a given position u along the 5th dimension), and the dilaton plays the same role as the onium 
— the corresponding wavefunctions are localized near the respective 'infrared' cutoff (u ~ uq 
for the dilaton and, respectively, r ^ Rq for the onium) and they exhibit power-law tails with 
exponent A in the 'ultraviolet' (large u and, respectively, small r). At this stage, it would be 
straightforward to repeat the analysis in Sect. 2 {mutatis mutandis) and thus determine the 
saturation momentum for the dilaton. However, before doing that, it is necessary to better 
understand the validity limits of the approximation in Eq. ( |3.9| ). We shall do that in Sect. 3.2, 
where we shall discover that, in the string theory context, Eq. ( p. 9] ) is less general than its analog 
in pQCD — it holds only within a limited range of values for r and p. These limitations have 
not been properly recognized in the previous literature. 

We conclude this subsection with a discussion of the interplay between the large- A''^ limit and 
the high-energy limit, which turns out to be quite subtle, for both weak and strong coupling. The 
tree-level string amplitude (corresponding to single Pomeron exchange) is of order oc 
In Eq. (|3.5D this factor has been canceled by a factor implicit in the 7^-current. (The 
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latter receives contributions from all the fields in the = 4 SYM theory which are charged 
under the 7?.-symmetry; these fields are scalars and Weyl fermions in the adjoint representation 
of the colour group, hence they have N"^ — 1 ~ N"^ degrees of freedom.) 

This is similar to what happens in QCD at weak coupling: the BFKL result ( 2.1C| ) involves 



an overall factor a1 ~ a'l/N'^ which is partially canceled in F2 by a factor Nc coming from 
the photon wavefunction, Eq. ( p.4D . Yet, the unitarity corrections apply directly to the dipole- 
dipole scattering (and not only to -F2), hence they become important when the small factor l/N"^ 
manifest in Eqs. ( p.lO[ ) or ( |2.15| ) gets compensated by the energy evolution. As discussed below 
Eq. ( p. 15 ), for the BFKL Pomeron this happens at a critical rapidity ~ {l/as)\i^{N'^/a1). 



For r ^ Tcr, multiple Pomeron exchanges become important and ensure unitarization. 

Returning to strong coupling and the dual string theory, there the multiple Pomeron ex- 
changes correspond to higher genus string amplitudes which are suppressed by powers of 
Hence, unitarity corrections become important at ~ InA'^^. We thus see that, in the strict 
large- A^c limit {Nc — > 00 at fixed energy), there is no issue of unitarity: higher loop ampli- 
tudes, being proportional to inverse powers of can be made arbitrarily small even at very 
high energy. On the other hand, when Nc is large but finite, the single Pomeron exchange 
has only a limited applicability, and beyond that one has to sum over an infinite number of 
multi-loop diagrams corresponding to arbitrarily many Pomerons exchanged in the t-channel. 
What is however different from the situation at weak coupling is that, for generic values of r 
and p ~ InQ^, the amplitude in the vicinity of the unitarity line is not dominated by the single 
Pomeron exchange anymore. This is why, in general, the saturation line cannot be inferred from 
the previous formulae in this section, but requires some further analysis. 

It will be our goal in what follows to determine the onset of the various types of unitarity 
corrections in the parameter space of r and p. To have access to the full structure of the high- 
energy 'phase diagram', we shall assume that Nc > e^ — so that diffusion effects (which, we 
recall, require a rapidity evolution r ^ \/A) become important before unitarization. But the 
opposite case e^ <C N^ will be included too in our results, as a special limit. 

3.2 Limit of the single Pomeron approximation : A lesson from a sum rule 

The fundamental difference between the physics of unitarization at weak and, respectively, strong 
coupling comes from the fact that, in the strongly-coupled gauge theories with a gravity dual, 
the high-energy scattering amplitude is predominantly real (it starts as one graviton exchange), 
while at weak coupling it is predominantly imaginary (in both QCD and AA = 4 SYM theory). 
Of course, the large real part at strong coupling is irrelevant for the total cross section so long 
as the latter is computed at tree-level. But at higher loop level, say, starting with two-graviton 
exchange, this gives rise to an imaginary part via diffractive processes — i.e., processes where 
the final states (string excitations or higher Kaluza-Klein modes) are separated by a rapidity 
gap. On the other hand, the imaginary part of the tree-level amplitude ( p.5| ) is associated 
with the inelastic production of excited string states, with no rapidity gap^. It turns out that, 
when p ~ InQ^ is larger than a certain critical value, the diffractive processes dominate over 



^Strictly speaking, the DIS process is always inelastic, since the virtual photon cannot appear in the final 
state. Here, we use the word 'inelastic' to suggestively distinguish the final states without rapidity gap from the 
diffractive ones, which exhibit such a gap. 
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the inelastic ones in the approach towards unitarity. We shall discuss this issue in detail in 
the next section, but for the time being let us proceed naively and study the onset of unitarity 
corrections on the basis of the tree-level results in the previous subsection, which include the 
inelastic processes alone. This will drive us into a paradox which will point out towards the 
necessity to include the additional, diffractive, contributions. 

After also including the essential factor the forward scattering amplitude ( |3.9D for 

the T^-boson-dilaton can be more suggestively rewritten as (compare to Eq. (|2.5| )) 

/•oo 

T{p,T) = / dp'To{p,p',T) mp')\^ (3.11) 
Jo 

where |^'(p')P = e"^''' and 

represents the elementary scattering amplitude between states localized around p and p' in the 
single- Pomeron-exchange approximation (the analog of the dipole-dipole amplitude in Sect. 2). 
The unitarity bound on this elementary amplitude determines the saturation line for a dilaton 
state localized at p'; namely, the condition Tq{p,p',t) ~ 1 for p = ps{T,p') implies 



PsiT,p') = p' + psoir), Psoir) = V^DriiOoT-lnN^), (3.13) 

where Pso{t) refers to a dilaton state localized at u ~ uq. (As in Sect. 2, we neglect the effects 
of the slowly varying prefactor 1/{2VttDt) in Eq. (p^.) 

Eq. ( 3.13 ) applies for t > Tcr = {l/u}o)lnN^ and should be compared to the corresponding 



result, Eq. ( ^.23 ), in pQCD: the similarity between these two results would become even closer 



if it was possible to extrapolate Eq. ( 3.13 ) to large rapidities r Tcr, where InA*"^ could be 
expanded out from the square root. However, as we shall later discover, such an extrapolation 
would be incorrect: in reality, Eq. ( |3.13| ) applies only within a limited range in t above T^r; 



namely, within the window < t < '^c; with — Tcr — 

2lnN^/VX (see Eq. (Q) below) 



This limitation is in turn related to a failure of our formula ( |3.12| ) for the elementary 
amplitude, which becomes incomplete at high energies and for large values of p. To illustrate 
this failure, we shall show that Eq. ( p.l2| ) leads to an unphysical result when used to evaluate 
the following moment of the structure function 

M{Q^)= / dxF2(x,Q2)= / dre-"F2(r,Q2). (3.14) 
Jo Jo 

By energy-momentum conservation, this moment should be independent of Q^] however, as 
we shall demonstrate in what follows, the (properly unitarized) amplitude for single Pomeron 
exchange, Eq. ( p. 12 ), fails to saturate the sum-rule at large values of Q^. 



Within the present approximations, the structure function F2 of the dilaton is given by 
Eq. ( |3.5D . To avoid unnecessary complications with the 7^-boson wavefunction, we shall here 
restrict ourselves to the structure function that would be measured by a single 7^-boson state 
localized at u (the analog of a single dipole projectile with fixed size r in pQCD). The coordinate 
u fixes the resolution scale, u"^ oc (cf. Eq. (|3.6[)), hence i*2(r, Q^) can be evaluated as 



F2{t,Q^) ~ iV^e^r(p,r) (3.15) 
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(up to an irrelevant prefactor), with p = h\{v?' /u^). The factor of A'^^ accounts for the color 
degrees of freedom of the fields contributing to the 7?.-current and e^ accounts for the factor 
in the relation (2.2) between F2 and the 'dipole' (here, 7^-boson) cross-section. 



We start by evaluating T(p, r) in the single- Pomeron approximation. This is given by 
Eq. (3.11) with the elementary amplitude from Eq. ( |3.12| ) corrected for unitarity violations. 



That is, for r > Tcr, we shall write (compare to Eq. ( |2.2'j| ) in pQCD) 

p—PsO 00 

T(/9,r) ~ j dp'^ e'^«"e-^Ti?^e-^^' + j dp'e"^^'. (3.16) 

As in Sect. 2, this representation is valid only in the regime where the amplitude is weak, 
T(p, r) <C 1, but it can be used to approach the unitarity limit T ~ 1, and thus compute the 
dilaton saturation momentum. The first integral in the r.h.s. of Eq. ( |3.16| ) involves the dilaton 
components living at relatively small values of u, which scatter only weakly off the incoming 
7^-boson. The second integral refers to the components at larger values of u, which undergo 
strong scattering, but are exponentially suppressed by the dilaton wavefunction. (Of course, 
when r < t^x only the first integral would be present.) Accordingly, the overall amplitude is 
dominated by the components living near the infrared cutoff, so like in pQCD. 

Indeed, the integrand in the first integral is peaked at p' = p — 2lS.DT] hence, for p < 2ADr, 
that integral is dominated by its lower limit /)' = 0, which yields 

T{p,t) ~ 4? e'^°^e-3i^ + 1 e-^^^-^^o) . (3.17) 

A 

This amplitude becomes of 0(1) when p ~ Pso{t)j meaning that the saturation momentum of 
the dilaton coincides with that of its component located at n ~ «o : /Os(t) = Pso{t)- However, 
the above calculation of ps is consistent only so long as ps < 2ADt, meaning psQ < 2ADt, 
which in turn requires t < Td, with 

In N'^ In N'^ 

rd = = ^7-^ • (3.18) 

u;o-^A2 1- (2/^/A)(A2 + l) 

When T < Td, one can check that the second term in Eq. ( |3.17] ) is exponentially suppressed next 
to the first one, and therefore 

1 p2 

T(/9,r) ~ — e'^o^e"4i>f when r < and p < 2ADt . (3.19) 

Thus, in this regime, the scattering is controlled by the component of the dilaton living near 
the IR cutoff uq = KE?, in full analogy to what happens in pQCD (recall the discussion after 
Eq. (|2.26D ). In particular, within the restricted interval r^r < t < r^, Eq. ( p. 191) can be rewritten 
in a form which resembles the expansion of the QCD amplitude near the saturation line (compare 



to Eq. (2.27)) 



T{p,t) ~ Q-C^-^^iP-P^) Q-'-^JB^ when < r < and ps < p < 2ADt , 

(3.20) 
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with the 'anomalous dimension' 7(r) given by 



1 - 7(r) 



Pso(t) 
2Dt 



uoQT-\nN^ 



(3.21) 



Note however that, even for p close to ps-, where the diffusion term can be ignored in Eq. ( |3.20| ), 
this amplitude does not show geometric scaling (unlike in pQCD), because the r-dependence of 
the 'anomalous dimension' ( 3.21| ) cannot be neglected within the corresponding validity range. 
This difference can be understood as follows: in pQCD, expressions like Eq. ( 2.27D are valid 
up to arbitrarily large values of r, much larger than the critical value ~ (l/os) ln(l/a^) for 
the onset of saturation; then properties like geometric scaling emerge via asymptotic expansions 



valid at T ^ By contrast, in the present context, Eq. ( 3.20 ) is valid only within a limited 
range of values for r above Tcr, where such asymptotic expansions are not possible anymore. 

The case where p > 2ADt can be similarly treated. For r < r^, the saddle point p' = 
p — 2ADt lies within the integration range and gives the dominant contribution to the overall 
amplitude: 



1 



,(c.o+DA2)r-Ap ^j^g^ ^^^^ P>2ADt, 



(3.22) 



in agreement with a corresponding result in Ref. ||T^. For r > t^, Eq. (3.16) would predict 
T{p,t) ~ e~^(P~P=o) ^ ^]^jg result cannot be trusted anymore since, as we shall shortly argue, 
Eq. ( |3.16| ) is not a reliable approximation for such large values of p and r. 



To see that, consider the prediction of Eq. ( |3.16D for the sum-rule (3.14). When computing 
the integral there, the amplitude ( |3.16| ) can be replaced by its simpler expression in Eq. ( |3.19| ) 
(this will be checked later). To also account for the unitarity bound, we shall write 



T{p,r) 



N? 



e 4Dt for r < Ts{p) 



for r > Ts{p) 



(3.23) 



Here, Ts{p) is the dilaton saturation line expressed as a function of p (the 'inverse' of the function 
Ps{t)); that is, Ts{p) is the rapidity at which the unitarity corrections become important on the 
resolution scale fixed by p. Using Eq. ( 3.13| ), one immediately finds 



InTVi 




+ 



(3.24) 



After using ( 3.23| ) together with the estimate ( 3.15| ) for F2, the integral in Eq. ( 3.14| ) becomes 



M{p) 



dre"(^"'^o)^e"4S57 + qP / dre~^ = Mi + M2 . 



(3.25) 



The second integral is trivially performed, and its result is found to rapidly vanish at large p: 



M2{p) ~ N^ePe'^^'-P^ — > exp 



AV4 
' V8 



p> for p —>■ 00. 



(3.26) 
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Hence, the only way to satisfy the sum-rule would be via the first integral, Mi. The correspond- 
ing integrand is strongly peaked near the saddle point at 



r*{p) = ^p. (3.27) 

Note that, for r ~ r*, we have 2ADt ~ Ap > p, which justifies our use of the simplified 
expression ( p.l9| ) for the single-Pomeron amplitude in Eq. ( |3.16| ). 

So long as T*{p) < Ts{p), this saddle point lies inside the integration domain and can be 
used to estimate the integral. The ensuing result for Mi{p) is independent of p, showing that the 
sum-rule is saturated by the single-Pomeron exchange, as it should for this to be a consistent 
approximation. This is the situation considered in Ref. jl^, within the framework of the strict 
large-A'^c limit. Note, however, that in that limit, the saturation rapidity was effectively 
pushed up to infinity and unitarity has never been an issue^. 

However, when Nc is kept finite (although large), unitarity is an issue, which modifies, in 
particular, the previous calculation of Mi{p) at large values of p. Namely, the function T*{p) 
grows faster than Ts{p), so the two curves cross each other at a point (see also Fig. |3|) 



Hence, when p > pc, the saddle point ( 3.27| ) lies outside the integration range. Then, the integral 



giving Ml is dominated by its upper cutoff at r = Ts{p), which however yields only a tiny result, 
of the same order as M2 (since for r ~ Ts{p), the amplitude ( 3.23| ) is 0{1)). 



Thus, for p > pc the sum-rule is not saturated by the single Pomeron exchange anymore, 
meaning that some new physical mechanism, which was missed by this approximation, should 
become important in this region. The above calculation gives us some indications in that respect: 
We have already noticed that, when increasing p towards pc, the saddle point ( |3.27 ) approaches 



the saturation line ( p.24|) , that it crosses exactly aX p = pc- It is then reasonable to assume 
that, for p > Pc, the saddle point should remain along the saturation line: T*{p) = Ts{p). (One 
cannot have T*{p) > Ts{p), since the integral over r in the saturation region is always dominated 
by its lower limit Ts{p), as manifest on Eq. ( p. 251 ).) If so, then the sum-rule at large p can 
be estimated as in Eq. (|3.26| ): M{p) ~ exp{/9 — Ts{p)}. This estimate, together with the 
condition M{p) = 0(1), allows us to predict the form of the saturation line at large p > pc- 

Ts{p) ~ p + InN^ , or Ps{t) ~ r - IniV^^ . (3.29) 

This discussion suggests that the complete saturation line should actually involve two branches, 
given by Eq. ( |3.24|) for p < pc and, respectively, Eq. (|3.29| ) for p > pc- Remarkably, these two 
branches appear to continuously match with each other at the transition point ( |3.28| ) : both 
their values and their first derivatives coincide at this point (see Fig. ^). 

What should be the unitarization mechanism responsible for the saturation line in Eq. ( |3.29D ? 
Once again, the previous sum-rule calculation gives us an indication in that sense: Along the 

^This is also the situation in pQCD, where the sum-rule (3.14) is dominated by relatively large values of 
X ~ 0(1), hence its calculation is insensitive to unitarity corrections, or even to the BFKL evolution. 
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line which controls the sum-rule, the energy dependence of the 7^-boson-dilaton amplitude is 
precisely that of the massless graviton {not Pomeron) exchange: 



= e(i-4/^)" = e" e-f . (3.30) 



(We have repeatedly used r = {^/X/A)p, cf. Eq. ( 3.27| ).) Eq. ( 3.30I ) is indeed the same as the 



propagator of a massless graviton. This result is quite suggestive : The shift of the Pomeron 
intercept 1 — > in Eq. ( ^ ) can be understood as a result of gravitons acquiring a mass due 
to the curvature of AdS^ (see Sect. 4.2 below for details). However, along the line Eq. ( |3.27| ), 
and in fact in the whole region on the right of this line as we shall see, the presence of massless 
gravitons is inevitable. The graviton exchange is dominantly real and much more non-local than 
the Pomeron exchange, which is non-local only due to diffusion. However, double, or multiple, 
graviton exchanges can generate imaginary parts, via the diffractive final states alluded to at 
the beginning of this subsection. As we shall explicitly see later on, such multiple exchanges are 



indeed responsible for unitarization at p > pc and they produce the saturation line in Eq. ( |3.29| ). 
4. Mapping the high— energy 'phase diagram' 

Motivated by the above considerations, in this section we shall study the Pomeron propagator 
in the AdS^ geometry, with due attention to both its real and its imaginary part. This will allow 
us to follow the transition between a genuine Pomeron behaviour at relatively low values of p, 
where the curvature effects are important and give rise to an imaginary part and to diffusion, 
and a massless graviton propagator at large p, which is long-ranged and predominantly real. 
Then we shall explain how multiple graviton exchanges generate an imaginary part via diffractive 
processes, and compute the corresponding saturation line. As a preliminary, we shall make some 
comments on the anomalous dimension of twist-two operators at strong coupling, which justify 
the use of the 'diffusion approximation' in the approach towards the saturation line. 

4.1 Anomalous dimension of twist two operators in AA = 4 SYM 

The variable 7 introduced in Sect. 2 has the interpretation as the anomalous dimension of the 
twist-two, spin-j gluonic operator 

Tr[F+(I)+)^-2F+'^], (4.1) 
with j analytically continued to non-integer values. (In the BFKL context one has j = 1 + x(7)j 



see 



3711 for discussions on the anomalous dimension and the operator product expansion in that 
framework.) Recall that the total dimension A(j) of the operator Eq. ( [4.lD is given by'^ 

A(j) = i + 2 - 27(j) . (4.2) 

In the weak coupling expansion, 7 ~ 0{as), while the BFKL resummation gives a number of 
0(1) (and comprised between and 1/2; cf. Sect. 2). On the other hand, for AA = 4 SYM at 



*This is the common definition of 7 in tlie smali-a; literature, wfiidi differs by a factor —2 from tfie one wfiicfi 
is most often used in tfie literature. A yet different convention 7 1 — 7 is sometimes used. 
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strong 't Hooft coupling A = g^N^ ^ 1, the AdS/CFT correspondence predicts that the twist- 
two operators with j ^ 2 receive a large, negative (in our conventions), anomalous dimension. 
Specifically, Ref. |Q has found 



for j » ^/\ and, respectively. 



7(j)«^-y|A3, (4.4) 



for 1 ^ j ^ \/A. These results were obtained by computing the energy and the angular 
momentum of a 'folded' closed string state rotating in the AdS^ space. When j <^ \/A, the string 
is spinning in a small region compared with the radius R, while the case j ^ \/A corresponds 
to a macroscopic string stretched over a distance larger than R. 

Because of the large anomalous dimension for twist-two operators (excepting the energy 
momentum tensor which has zero anomalous dimension 7(2) = 0), it has been argued in [^] 
that the OPE in DIS at strong coupling must be considerably reorganized. The dominant 
contribution comes from the 'protected' (i.e., non-renormalized), double-trace, operators which 
create and annihilate the entire hadron of interest. These operators are nominally higher-twist, 
but since their dimension is of order unity, their contribution is less suppressed compared to the 
twist-two operators having \ ^ 1 for generic j 7^ 2. Hence DIS at moderate values of x 
can be pictured as a process where a virtual photon knocks off an entire charged hadron - a 
drastically different picture from pQCD where a photon knocks off partons inside the hadron. 

While plausible in the Bjorken limit (Q^ ^ 00 at fixed r), this is no longer the case in the 
Regge limit (r ^ 00 at fixed Q^). The important regime at high energy is j ~ 2. Since 7(2) = 
exactly, 'y(j) remains small in the vicinity of j = 2, and the twist-two operators (analytically 
continued to complex values of j near j = 2) should still give the leading contribution. Notably, 



Refs. |38] and |]1^] arrived at the same expression for in the 'diffusion approximation' from 



very different arguments. It is obtained as the solution to 



where 



JO = 2 - = 1 + , (4.6) 

is the Pomeron intercept. Notice that the property 7(2) =0 is indeed verified by the solution 
to Eq. ( f4.5[ ). Some light on the origin of Eq. ( [1.5| ) will be shed by the manipulations in the next 
subsection. The inverse function j{'y) is approximately quadratic in 7 and has a minimum at 

70 = 1-^. (4.7) 

with the minimum value ^(70) = Jo- 
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In general, the diffusion approximation is expected to be valid only in the vicinity of this 
minimum, i.e., when 7 ~ 70, j ~ jo- Away from j = 2, one could a priori expect significant 
nonlinear corrections to the right hand side of Eq. ( [4.51 ). This is precisely what happens at weak 
coupling, in the context of the BFKL equation, where the diffusion approximation corresponding 



to 7o = 1/2, cf. Eq. ( 2.14 ), significantly deviates from the exact BFKL characteristic function 
already when 7 = 7^ ~ 0.37. (This is why, in that context, it would be incorrect to compute the 
saturation line by 'unitarizing' the amplitude ( ^.15 ) obtained via the diffusion approximation). 



However, it turns out that at strong coupling the diffusion approximation is much more 
robust and has a much wider validity range than at weak coupling. For instance, one can see that 



Eq. (4.5) correctly reproduces the leading behavior, Eq. (4.4), for large (but not asymptotically 



large) j. (As a matter of fact, Ref. |35] has determined the proportionality constant in Eq. ( |4.5| ) 
in such a way to match Eq. (4.4). But the derivation of Ref. [jl6| does not rely on this matching.) 
Moreover, Eq. ( [4.5[ ) coincides with the well-known anomalous dimension formula of a scalar field 
in AdS^ 

A = 2 + V4 + m2i?2 , (4.8) 

where m is the mass of the state lying on the leading Regge trajectory (i.e., the highest spin 
state in a given excitation level n = 0, 1, 2, . . . ) : 

2 4n 2(j-2) ,^ , 

o! a 

This shows that, in this strong-coupling context, the quadratic relation j ~ 7^ ('diffusion 
approximation') has its roots in the Regge behavior of string excited states ~ m? ~ j, which 
is expected to hold in AdS^ for generic values of j up to j ~ ^/A (see Eq. This range is 

large enough to cover all values of j and 7 of interest at high energy, as we shall later verify. 

4.2 The Pomeron propagator 

In Eq. ( |3.5| ), the Virasoro-Shapiro formula in flat space was used rather heuristically ||l^. For 
a more coherent discussion on both the real and the imaginary parts of the amplitude, it is 
important to carefully derive this amplitude from the fundamental S'-matrix in string theory, 
starting with the operator representation of the latter in curved spacetime. Up to a factor g1 to 



be introduced in the final result, this representation reads [p3| , 1£] 



d^z (V(oo)V(l)W(z)>V(0)) f« / d^z (v(oo)V(l)z^""^z^""^>V(l)>V(0)\ 

J\z\<i \ ' 

5r 



/v(cx))V(l) Y° W(1)W(0)^ • (4-10) 
\ Lo + Lo-2 I 



Here, V 's and W 's are the vertex operators corresponding to dilatons and, respectively, gauge 
bosons in the (0, 0) picture, and the integration is over the string world-sheet with spherical 
topology. The average is understood in the sense of the string path integral, which involves an 
integral over the zero modes and a path integral over the nonzero ones. Furthermore, Lg and 
-Lq are the right- and left-moving Virasoro operators which generate the scale transformations 
on the world-sheet. The kernel £^J{Lq + Lq — 2) in the last expression can be recognized as 
the propagator of the closed string exchanged in the t-channel. 
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One can diagonalize the operators Lq and Lq by inserting a complete set of string states — 
the 'Pomeron states' — , which are the states exchanged in the t-channel at high-energy [|l^]. 
To that aim, one inserts the identity as a sum over the respective projection operators 

where the sum runs over all the positive, even, values of the spin : j = 2, 4, corresponding to 



the string states on the leading Regge trajectory in Eq. (4.9). Alternatively, the sum over j can 
be realized as an integration in the complex j-plane, as indicated in the r.h.s.; the contour C 
encircles the positive part of the real axis and is conveniently deformed into a contour L which 
runs parallel to the imaginary axis, on the left of all the poles at j > (see Fig. ^ below). (More 
precisely, this contour deformation within the Sommerfeld-Watson representation is truly an 
analytic continuation from the t-channel amplitude related to the physical, s-channel, amplitude 
via crossing symmetry {t ^ s); see, e.g., ||40|| .) 

In the equation above, is the 'Pomeron vertex operator' which is dual to the twist- 
two operator, Eq. ( [4.1[ ), of the gauge theory (see also [0). For large u, it takes the form 

~ W-^'^'''{d,X^d,X^)-2 . (4.12) 

Capital letters denote world-sheet fields, while lower case letters refer to their zero modes. Also, 
light -coordinates are defined in the standard way: = {X^±X^)/V2, with the collision axis 
taken along X^ . For what follows it is useful to notice that the vertex operators associated with 



the dilaton (V) and the 7^-boson (W) admit representations similar to that in Eq. (4.12), but 
with the [/-factor replaced by the corresponding wavefunction — ^(U) in the case of the dilaton 
and, respectively, AmiU) for the 7^-boson. 

The operator (|4.12| ) has total dimension A = 2 + 2iu; together with Eq. ( [4. 2]) , this implies 
the following anomalous dimension : 

7 = |-ii/, (4.13) 

which appears to be different from the one usually associated with a leading-twist operator 
in the diffusion approximation, cf. Eq. ( p.l3| ). Let us open here a parenthesis to explain this 
discrepancy: in the leading logarithmic approximation used throughout Sect. 2, the reference 
scale of the logarithm r = Ins is ambiguous and usually taken to be a product of the two 
external scales. This ambiguity is resolved in the next-to-leading log approximation (NLLA) 
where 7 is shifted in the following way [0] 



\QAJ J XQ"^) 



(4.14) 



to ensure that the Pomeron behaviour {s/Q'^y^^ matches with the definition of the Bjorken-x 
variable in DIS (at high energy). 

The action of Lq on the Pomeron vertex operator has been evaluated for j ~ 2 |l^ as 

ioVj, = Uvt. = (i - f A|) V*. = (I + ^i±i) Vt. . (4.15) 



-25- 



where (c.f. Eq. 



(4.16) 



The z/"-integral (or the ^-integral) sets Lq = Lq = 1, which after also using Eqs. ( [4.13|) -( ^.15| ) is 
recognized as the previous condition ( |4.5[ ). So, Eq. ( [4.5| ) simply states that is a (1,1) tensor 
on the world-sheet in the curved background of AdS^ p^]. 

We are now prepared to evaluate the string amplitude in Eq. ( [4.10 ) in the semi-classical 
approximation. After using Eqs. (4.11) and ( 4.1^ ), one finds 

V(oo)V(l) — W(1)W(0)\ = 







/ 



dj du 1 + e" 



■ITT J 



1 



47ri 



sm^J j-jo + 



(V(oo)V(l)V^r_J(V+,>V(l)W(0)) . (4.17) 



When evaluating expectation values in the semi-classical limit, one can replace an operator 
F{U) on the world-sheet with its value F{u) at position u (the zero mode U). Besides, the 
path integral over U reduces to the ordinary integral over u, with the appropriate, invariant. 



measure. Since the string propagator in Eq. ( 4.17 ) is non-local in u, some care is needed when 
constructing the measure. Namely, in the decomposition of unity, Eq. ( 4. Ill ), we need to insert 
the vertex operators at different positions u and u' . The appropriate decomposition reads then 

djl + e"''^^' r dv f du' 
2^ 



1 



^ 2i smnj 

where G is the determinant of the AdS^ metric Gmn, and is the respective metric compo- 
nent; Eq. (3J.) imphes G~^~ = /u^ and ^/—G = (u/R)^. Eq. (|4.18| ) follows from the following 
decomposition of the (5-function, which can be easily checked: 



R 



V^{G' 



+-V 



V,--.(t^'))(V+,(n)| 



(4.18) 



dn' 5{u 



R ^ ' 2t:\R 



r) 



(4.19) 



One can now insert the explicit expressions for the vertex operators in the u-representation and 
then perform the integral over v by deforming the contour in the complex plane. Using the fact 
that u > u' in DIS at large (recall that u Q, whereas u' ~ uq), we shall close the contour 
in the lower half plane, and thus pick up the pole at 



This gives 

duV^ 



dj 1 + e-'^^ (a'^s^y 



(j - Jo)- 



(4.20) 



2i sinvri VJ - Jo 

X / dn'^^(G"+-)^'$t$(^')^'i-2+2iKi)^i-2-2Mi)^^(^^^)^m^^^ 



+ 



dj 1 + e" 



ITT J 



2i sinvrj Vj - Jo 



du 



u 



^<i>t$(n') 

u 



/2\ l-7(i) 



A^{u,q)A"'{u,-q) + --- 



(4.21) 
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where the dots stay for terms which are subleading at high energy. The factor (cf. Eq. ( |3.6| )) 
arises after contraction with vertex operators which carry large momenta [16|. Apart from the 



external wavefunctions, Eq. ( 4.21D is symmetric under the exchange of u and u' as it should. 



After including the factor gl ~ l/N"^ and substituting the solution of the Maxwell equation 
for Am{u) one can identify the forward scattering amplitude for the current-dilaton scat- 
tering, whose imaginary part yields the F2 structure function: F2 oc N^Q^lmT. [We use the 
notation T for the complete scattering amplitude, including both the real and the imaginary 
part; hence, T = ImT.] One finds 

\6 



^<^,i^u')^^ 4 . (4.23) 

u' ^ ' J 2i sinvrj ^rj^^\y2 ) ^ ' 



Eq. ( |4.23 ) exhibits 'radial factorization': it is expressed as a convolution in u of pieces which 



can be recognized as wavefunctions for the 7?.-boson and, respectively, the dilaton, and a kernel 
which is by construction the Pomeron propagator in string theory — a sum over all string states 
exchanged in the t-channel at high energy — , and which after multiplication by gl plays the role 
of a scattering amplitude between states localized at u and u' . Clearly, this u-factorization is 
the string counterpart of the 'fc^-factorization' valid in the high-energy regime at weak coupling 
(as exhibited, e.g., in Eqs. ( |2.3| ) and (^^)). 

Let us concentrate on the elementary Pomeron amplitude, that we would like to compare 
to its previous approximation in Eq. ( |3.12| ). With our usual notation p = \i\{u'^ / Uq) , this can be 
rewritten as 



Tp(/9, p',t) = — — — : : 7== exp { {p - p') 

Jl 2i smvrj VJ - Jo 




(4.24) 



where we have used Eqs. ( [4.13 ) and ( [4.20 ) to explicit the function The contour L runs 



parallel to the imaginary j axis and crosses the real axis between the branch point at j = jo ^ 
2 — 2/\/A and the pole at j = 2 (see Fig. ^. 

For sufficiently high energies, the integral can be evaluated by deforming the contour L L' , 
with L' encircling the branch cut that ends at jo, as shown in Fig. ^. Alternatively, and more 
conveniently for our present purposes, one can use the saddle point method. Namely, the 
integrand has a saddle point at 

/-.o + ^%^. (4.25) 

where we have used r ^> p — p', as appropriate at high energy. So long as \/A(p — p')/4 < this 
saddle point lies between jo and 2, and then the use of the saddle point approximation is indeed 
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Figure 2: Tiie analytic structure of the integrand in Eq. (4.24) in the complex j plane, together with 
various contours used when evaluating the integral. 



justified. By also using e" one thus finds that the high-energy amphtude 

develops an imaginary part^ Tp = ImTp, with 

Tp{p, p', r) « 4? e^°" e-^^-^°'>^f-f"> e~^^^ , (4.26) 

where = 1 — 2/ a/A, 70 = 1 — 1 / "/A, and D = 2/ ^/X. 

Eq. ( 4.26| ) coincides with the formula previously used in Sect. 3, cf. Eqs. ( ^^ ) or ( 3.12| ), up 
to a factor e~^^~^°^^ = (u'"^ /u^^^^ . This additional factor has the same origin as the factor 

case, Eq. ( ^.15| ) — namely, is reflects the high-energy 'anomalous 
dimension' — , but it is less important here because the anomalous dimension 70 is close to one 
(in contrast to the BFKL case where 70 = 1/2). And indeed, by inspection of Eq. ( 4.26 ), one 
can check that the additional term (l—jo){p — p') in the exponent remains negligible so long as 
p — p' <^ and hence it cannot affect our previous conclusions in Sect. 3.2. Note also that, in 
this high-energy approximation, the graviton pole ~ 1/t is absent in the real part of the string 
amplitude — in contrast to what happens in flat space — , because of the shift in the Pomeron 
intercept by the space-time curvature: sinvrjo / 0. 

However, the graviton pole reappears when moving to sufficiently large values of p (for a 
given r), as we explain now. Indeed, when 

^{p-p')>T, (4.27) 

the saddle point j* exceeds 2, and then one has to deform the contour from C to C' (see Fig. 
0), in order to separately evaluate the pole of (l/sinvrj) at j = 2. This leads to 

THp,p'.r) . I.e-I-"') + -1 / P±£:i^^i^-.m^-^^ . (4.28) 

J 2i smvrj VJ - Jo 
c 

^The Virasoro-Shapiro formula in flat space also gives an imaginary part in a similar manner. e~"^" ^^'^ /t ~ 
1/t — ina' /2. This is how Eq. (3.5) was actually derived. 
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The first, purely real, term tx e^{u''^ /u^) is recognized as the propagator of the elementary, 
massless, graviton in AdS^. For large r and p (such that the condition ( f4.27| ) is satisfied), 
this pole term dominates over the remaining contour integral. Indeed, the function /(j) = 
{j — + ^{j)[p — p') obeys /(2) = r, f'{j*) = 0, and f'{j) < for any real j within the interval 
2 < j < j* ■ Hence, whenever ( [4.27| ) is satisfied, one can choose the contour C (see Fig. |2|) to 
cross the real axis at some ji with 2 < j'l < j* , and then Re /(j) < r for any j along C. 

Thus for high energy and sufficiently large p, the 'Pomeron' exchange is dominated by the 
elementary graviton and the amplitude is predominantly real. Since the graviton propagator is 
very non-local in u, it can connect the large distance between u ~ QR^ and u' ~ AR^ without 
much suppression. Accordingly, with increasing r, the corresponding amplitude becomes of 0(1) 
relatively fast (much faster than for the diffusive Pomeron in Eq. ( [4.26 ) !), namely when 



p-p' + lniVi, (4.29) 



in agreement with Eq. ( 3.29| ). Of course, being real, this single graviton exchange does not 



contribute to F2. Yet, multiple such exchanges can do so, as we shall shortly explain. 

To summarize, the single-Pomeron-exchange approximation to F2{x,Q'^) is expected to be 
a legitimate approximation at high energy and strong coupling so long as /) = In Q^/A^ < 4r/\/A 
and for rapidities r < Ts{p), with the saturation rapidity Ts{p) given by Eq. ( p^ . The two 



delimitating curves intersect with each other at the point {pc,Tc), Eq. ( 3.28| ). This is also the 
point where starts the large-p branch of the saturation line, cf. Eqs. ( |3.29[ ) or (4.29), to be 
further discussed in the next subsection. 

4.3 Saturation line from multiple graviton exchanges 

The previous analysis shows that the physics of unitarization in DIS at strong coupling and 
relatively large is conceptually different from what happens either at weak coupling, or at 
strong coupling but for lower values of Q^. In the last two cases, which are the more standard 
ones, the black disk limit for the scattering amplitude is approached by unitarizing the single 
Pomeron exchange, which in turn requires energies which are high enough for the two scales 
involved in the collision to be joint by diffusion. But for strong coupling and large Q^, the single 



Pomeron exchange is unable to saturate the energy-momentum sum rule ( 3.14 ) — essentially, 
because its imaginary part decays too fast at large Q^. However, precisely in that large-Q^ 
region, the full Pomeron propagator in string theory develops an additional, real, part, which 
describes elementary graviton exchanges and has the potential to generate stronger interactions, 
since very non-local in u. This new component can contribute to the DIS structure functions too, 
but only via multiple scattering. In this subsection, we would like to argue (without performing 
a fully fledged calculation of multiple scattering in AdS^) that multiple graviton exchanges 
dominate indeed at high energy and large (in the vicinity of the saturation line) and, in 
particular, they saturate the sum rule in the way already suggested after Eq. ( 3.26]) . 



The presence of a large real part in the Pomeron propagator describing the massless graviton 
exchange is a common feature in string theories, and hence also in gauge theories with a gravity 
dual. Multiple scattering in such theories is described by multi-loop diagrams, and at each order 
of this loop expansion the dominant contribution is provided by the exchange of n gravitons 
(for a (n — l)-loop amplitude). It is well-known that in fiat space such higher loop terms 
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exponentiate in the eikonal approximation [|^, 12 1, and the same was recently shown to hold for 
high energy scattering in AdS^ |18, 17|. The eikonal approximation breaks down at very small 
impact parameters and/or at extremely high energies (see below). In a string dual description 
of DIS at strong coupling, the impact parameter in the fifth dimension is large u/u' ~ Q/A, and 
therefore the eikonal resummation is expected to work reasonably well. 

In general, the exponentiation is done in the three-dimensional impact parameter space 
spanned by the radial coordinate u and the two-dimensional, physical, impact parameter b (ne- 



glecting nonlocality on S^), whereas our previous formulae, like Eq. ( 4.24 ), are already integrated 
over b. While it is possible to undo this integral and analyze the three-dimensional propaga- 
tor following fl^, one may anticipate, following the example of QCD at weak coupling (cf. 



Sect. 2; see also Ref. |44]), that a simplification occurs for relatively central impact parameters, 
6 <C 1/A. DIS is characterized by the large resolution scale ^ A'^ which determines the size 
of the 'active' region in the transverse space. Hence, so long as 6 ^ 1/A one can neglect edge 
effects and work with scattering amplitudes averaged over b, so like Eq. ( |4.24 ). 

In the present context, the eikonal approximation amounts to exponentiating the tree- 



level 'Pomeron' amplitude Eq. (4.24), which thus becomes the phase shift in the 5'-matrix: 
S = expjiTp}. In the interesting domain at large p = ln(Q^/A^), namely, p > 4r/\/A, the 
Pomeron amplitude is given by Eq. ( 4.2g| ) and has both real and imaginary parts. 



The imaginary part describes 'inelastic processes' in which the final states are (generally 
massive) string excitations, as obtained by cutting through the Pomeron. This imaginary part 
can be evaluated in the saddle point approximation, with the result displayed in Eq. ( f4.26| ). 
By itself, this contribution would generate the saturation line in Eqs. ( 3.13| ) or ( p. 29)) , which 
however has been argued to be inconsistent at large p > pc, in Sect. 3.2. 

The relevant contribution at large p is rather the one generated by the real part of the phase- 
shift, that is, by the exponentiation of the elementary graviton exchange. The corresponding 
imaginary part starts at the level of the two-graviton exchange, where it describes, in particular, 
the cross-section for the elastic scattering between two states localized at u and u' , respectively. 
Note however that, in the context of DIS, this process is not truly elastic — rather, it is a 
particularly simple diffractive process, which contributes to F2 — , because the incoming virtual 
photon (or 7^-current) does not appear in the final state. The amplitude for such a double 
graviton exchange is estimated as 



T2g{p,p,T) 



p') 



(4.30) 



and becomes of 0(1) along the saturation curve (4.29). More general diffractive processes — 
characterized by a rapidity gap between the products of the collision — are obtained by 'cutting' 
the multiple scattering series in between successive graviton exchanges, and by allowing for more 
general diffractive states: the strings in the final states need not be exactly the same as the 
incoming ones, rather they can be internal excitations of the latter, like higher Kaluza-Klein 
modes. Whenever this happens, the process is genuinely absorptive, already at the level of the 
elementary p — p' amplitude. 

The diffractive excitations of the string have been first discussed in Ref. in fiat space. 
Leaving a rigorous treatment of this effect in AdS^ for future work, here we only give an estimate 
of its magnitude. The internal excitation of incoming strings is due to the finiteness of the string 



-30- 



size. This effect can be understood as the result of the tidal force [45| experienced by one of 
the strings traversing a gravitational shock wave created by the other. In the operator eikonal 
approach of ||l^ , the real part of the phase shift reads 



,'2 



e'^ u'^ I / e' u" 



(4.31) 



where again capital letters denote operators, decomposed between zero and non-zero modes as 
U = u + U . The non-zero mode U is typically of the order the string size, as measured by a 
local inertial observer: 



U 



I - 



(4.32) 



so the expansion parameter in Eq. (4.31) is ll/R? = 1/\/A. The strength of string excitations 
leading to absorption is measured by the second term of Eq. ( [4.31| ). It reaches unity when 



(4.33) 



which, once again, is parametrically the same as Eq. ( [4.29| ). (Remember that we are in the 
regime where InA'^^ > A ^> In A.) We thus see that all types of diffractive processes — both the 
quasi-elastic ones, cf. Eq. ([4.30| ), and the genuinely absorptive ones, cf. Eq. ( 4.31 ) — contribute 
on equal grounds to the imaginary part of the forward scattering amplitude, and approach the 



unitarity limit along the curve exhibited in Eq. ( 4.29 ). This is in contrast to the flat space result 
1 12 1 where the imaginary part due to string excitations is suppressed with respect to the real 
part by an inverse power of the impact parameter l1 /b"^ . 

So far, the diffractive unitarity line ( [4.29 ) has been derived for the scattering between two 
states localized at p and p' . But it is easy to check that a similar line, with p' = 0, applies to 
the dilaton as whole. Indeed, after performing the convolution with the dilaton wavefunction, 
the relevant amplitude can be evaluated as in Eq. (3.16), which gives 



P-PsO 



?2s(p, t) 



dp' 



dp e 



-Ap' 



(4.34) 



p-PsO 



where now pso = t — InNl, cf. Eq. (|I29| ), and we have assumed that p > pso- In writing 
the equation above, we have focused on the two-graviton exchange, cf. Eq. ( 4.3C| ), since this 
is the dominant contribution when approaching the saturation lien from large values of p. At 
this point, it is important to remember that the dilaton conformal dimension satisfies A > 2, 
as manifest on Eq. (|4.8|). Hence, the first integral in Eq. (4.34) is dominated by its lower limit 



p' = and, moreover, the second integral is exponentially suppressed w.r.t. the first one. Thus, 
once again, we find that the scattering is controlled by the dilaton component living near the 
infrared cutoff at n ~ tio- This yields 



T2g{p,T) 



1 2 



(4.35) 



which becomes of 0(1) for p ~ PsiT) = r — InN^, in agreement with Eq. (|3.29|). 
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For p > pc, with pc defined in Eq. (3.28), the curve ( 3.29| ) hes below the unitarity hne 
( 3.24| ) for inelastic processes (see also Fig. |3|), meaning that, when increasing r in this regime, 
the unitarization of DIS proceeds via diffractive processes. Therefore, Eq. ( |3.29D represents the 
physical saturation line for p > pc- This conclusion is further corroborated by the following facts: 
(i) at p = Pc, the two branches of the saturation line, Eq. (|3.24|) and, respectively, Eq. (|3.29|) , 
continuously match with each other, and (ii) for p > p^, the energy-momentum sum rule 
( 3.14 ) is saturated by the diffractive processes (see below), showing that no important physical 
mechanism has been left out by the previous analysis. 

To evaluate the sum-rule ( 3.14| ), notice first that a contribution to F2 due to the exchange 
of n gravitons (with n > 2), is parametrically 



F2(r,g^) ~ iV^T.,(p,r) ~ iV^Q! (^__^__^ . (4.36) 

When this is inserted within the Eq. ( 3.14| ), it is clear that, for any n > 2, the integrand is 
strongly peaked at the saturation value Tg = p + lnN^: indeed, the integrand rises exponentially 
with T so long as r < r^, but it is exponentially decreasing at r > r^. Moreover, at the peak 
value, Tng{p,T = Tg) ~ 0{1), so the integral is evaluated as 

dre-"F2(r,Q2) ^ |! g-^W ~ 0{l). (4.37) 



/ 



This result is independent of Q^, as anticipated. Interestingly, the above calculation also shows 
that, at high Q^, most of the total energy (or longitudinal momentum) of the hadron lies along 
the saturation line ( ^.29 ), that is, at relatively large rapidities, or small values of x. This is in 
sharp contrast with what happens in QCD at weak coupling, where the energy of a fast moving 
hadron is predominantly carried by the hadron constituents ('partons') with relatively large 
values of x, or small values of r, well below the corresponding saturation value Ts{p). We shall 
return to this observation in the final, discussion, section, where we shall attempt a partonic 
interpretation for our results at strong coupling. 

We conclude this section with some estimates for F2{t, Q"^) in the vicinity of the saturation 
line ( 3.29 ), to be compared to the corresponding results in pQCD, cf. Eqs. ( 2.3[1| ) and ( |2.31| ). 
Note first that Eq. ( p9D implies QHt) = h? {e" / N"^) . 

For relatively low virtualities, ^ Qg(r), one finds 



/ S V 

F2(r,Q2)^^ when Q'' < QI{t) , (4.38) 

which shows saturation, in the sense that F2 does not increases with r, and that F2/Q'^ is not 
divergent as — > 0. (A physical interpretation for this behaviour will be proposed in Sect. 5.) 

For much larger than Q1{t), the structure function is dominated by the two-graviton 
exchange and scales as 

^^(-«V^(|^)\ (4.39) 

or 

F2ir,Q^) _ NUQlir)^^ when » Q^r) . (4.40) 



A2 V Q 



- 32 - 



This is geometric scaling (compare to the weak-couphng result in Eq. ( |2.30| )), with an effective 
anomalous dimension 7^ = —1, which however is no longer associated with the anomalous 
dimension of the twist-two operator. Indeed, the contribution exhibited in Eqs. ( [4.39 ) or ( [4 .401 ) 



is obviously of higher twist order, and in fact this process has an analog in standard DIS: the 
'diffractive dissociation' of the virtual photon, 'y*p — > qqp, which involves the elastic scattering 



between the 'color dipole' (the qq pair) and the proton (see, e.g., [|^, 34, |4^). There are however 
some interesting differences between these results and the corresponding ones at weak coupling: 

(i) In pQCD, the diffractive scattering represents only a tiny fraction of the total DIS cross- 
section, whereas in the present, strong-coupling, case this has been found to be the dominant 
component at high energy, which in particular is responsible for unitarization. 

(ii) The diffractive structure function F<P in pQCD at high energy (small x) and large 

is formally a leading-twist quantity, /Q'^ ~ although the relevant dipole amplitude is, 

of course, of higher-twist order. This is so since, in computing , the convolution with the 



virtual photon wavefunction, Eq. (2.4), is dominated by the relatively large qq excitations having 



r » 1/Q (the so-called 'aligned-jet configurations') |34, This does not happen, however, 
in the present context at strong coupling, where the convolution in Eq. (4.23) is controlled 
by li ~ Uc ~ Q- Accordingly, the ensuing structure function is a higher-twist quantity, cf. 
Eq. ( 4.40| ), so like the elementary scattering amplitude in Eq. ( 4.30| ). 



5. Physical discussion: Towards a partonic picture at strong coupling 

Let us first summarize the main new results in this paper, before we propose a physical inter- 
pretation for them. (These results are also graphically summarized in Fig. ^.) By studying 
deep inelastic scattering at strong 't Hooft coupling and high energy within the dual string 
theory, we have inferred the position of the saturation line which separates the weak-scattering 
regime from the strong-scattering one in the kinematical plane t — p, with r = ln(s/(5^) and 
p = ln((5^/A^). This line is in fact the juxtaposition of two curves, matching with each other 



at the point {pc,Tc) defined in Eq. (3.28|), which correspond to the two different mechanisms 



responsible for unitarization in different kinematical regions : 

• For /3 < pc — 4 In / \/A, the scattering amplitude in the vicinity of the unitarity limit is 
controlled by a single Pomeron exchange; this amplitude becomes of 0(1) on the saturation 
line given by Eqs. (3.13) or ( |3.24| ), i.e., when = Q^(r), with the saturation momentum 



Q2(r) = A2e\/^^"('^«"-l'^^') for r^r < r < , (5.1) 

with LOO = 1- 2/\/A, D = 2/VX, and — In A'^ is the critical rapidity for the onset of 
unitarity corrections. 

For p > pc, the dominant contribution to the DIS structure function before unitarization 
comes from multiple exchanges of elementary gravitons; the reaches the unitarity limit on 
the saturation line shown in Eq. ( |3.24 ), i.e., for = Q1[t), with 



Q'W=^'^ for T>T,. (5.2) 
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Figure 3: Proposed 'phase diagram ' for DIS at high energy and strong couphng when In ^ VA- The 
continuous, thick, hne represents the saturation hne. The other boundaries are explained in the text. 



In Fig. ^, the lower branch of the saturation Une, Eq. (^), is also represented, as a dotted 
line, for r > Tc, where it indicates the 'unitarity limit' for inelastic processes, i.e., the line along 
which the contribution of the single-Pomeron exchange to the DIS amplitude becomes of 0(1). 
Note that the 'diffusion approximation' that was implicitly used in deriving the equation of this 
line in Sect. 3.2 is still valid^ for such large values of p because, along this line, j ^ 3, with 
the upper limit j = 3 corresponding to p — > 00. (This can be checked by using Eq. ( [4.25| ) 
together with the large-p version of Eq. (|3.24| ), that is, Ts{p) ~ {\^/'^p/\/2 + lnN^)/2.) Below 
that line and for p > pc, the DIS cross-section is dominated by the diffractive processes (whose 
contribution becomes of 0{1) already along the proper saturation curve, Eq. ( |5.2| )), whereas 
above that line, both diffractive and inelastic processes contribute on equal footing and together 



give an amplitude T = 1. Also, below the saturation line (5^), the DIS cross-section for p > pc is 
predominantly diffractive, of higher-twist order, and it exhibits geometric scaling, cf. Eq. ( |4.40| ). 
The upper corner of Fig. |3| denoted as 'strong gravity' will be discussed towards the end of this 
section. 



^As discussed in Sect. 4.1, the diffusion approximation is expected to be correct so long as j <^ y/X. 
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So far we have assumed N'^ > e^, to allow both regimes alluded to above to coexist with 
each other. But it should be clear that the opposite case <C e^ (with ^ 1, though) is 
included too in our results, as the special limit pc — ^ : then, the first regime squeezes to zero, 
meaning that there is not enough phase-space for the genuine Pomeron behaviour to appear — 
the amplitude reaches values of order one via multiple graviton exchanges already before the 
diffusion becomes important. 

Note that the above results do not involve special properties of the dilaton target, like its 
mass or its conformal dimension A. We therefore expect them to universally apply (at strong 
coupling) to all hadrons which, like the dilaton, have a dual state localized near the infrared 
cutoff no = Ai?^ in the radial dimension of AdS^ . The same universality applies to the form of the 
structure function F2{t, Q^) in the vicinity of the saturation line. In the weak-scattering regime 
at » Q1{t), F2 shows a rapid rise with r and a rapid decrease with Q^, F2 oc (e^'^/Q^), cf. 
Eq. ( 4.39| ), which is the hallmark of the double-graviton exchange. Furthermore, at low momenta 
(r), one finds F2 ~ N^Q'^/A'^, cf. Eq. (|]3|), a structure which is an almost automatic 
consequence of the definition ( 3.15| ) of F2 together with the black disk limit for the scattering 
amplitude. But in QCD at weak coupling at least, this structure is also tantamount of parton 
saturation, which makes it tempting to propose a similar interpretation in the strong-coupling 
case as well. 

Let us first briefly recall the situation in pQCD (see, e.g., the review papers [27| for more 
details): there, the structure function F2{t, Q^) measured in DIS can be interpreted as the quark 
distribution in the target, as measured in the special frame (the 'infinite momentum frame', 
or IMF) in which the target has a very large longitudinal momentum P ^ M. The 'quark 
distribution' is the number of quarks per unit rapidity which are localized in impact parameter 
space within an area ~ l/Q^ fixed by the resolution of the virtual photon; alternatively — via 
the uncertainty principle — , this is the number of quark excitations with transverse momenta 
k± ^ Q. Hence, by knowing F2{t, Q^), one can estimate the quark occupation number as 

where the derivative in the r.h.s. is evaluated at Q'^ = k'j_, Rq is the radius of the hadron disk, 
and the factor 1/Nc enters since we consider quarks of a given color. (Other numerical factors, 
like the number of spin states, or the number of flavors, are not explicitly shown.) Then, by using 
the estimates ( p. 30 ) and ( |2.31| ) for F2, one can distinguish between two physical regimes: (i) a 
dilute regime at high momenta k± » Qs{t), where nq{T,k±) scales like'^ Ug ~ {Q'i{T)/k'j_)^~'^'' , 
and thus rises rapidly with r due the BFKL evolution, and (ii) a saturation regime at momenta 
k± ^ Qs('7"), where Ug saturates at a value of order 1. This is quark saturation and, as explained 
in Sect. 2, it ultimately reflects the saturation of the gluon occupation numbers in pQCD at 
small X m, m, H, I, i : ng(r, A:^) ~ 1/a, for k^ < Qs{t). 

Returning to the strong-coupling problem of interest, it is clear that a straightforward par- 
tonic interpretation is now hindered by the lack of a simple relation between the operators which 
enter the structure of the 7^-current (for the M = 4 SYM theory, this current receives contri- 
butions from the Weyl fermions and from the scalar fields) and the physical excitations in the 

^For very large fc^i the anomalous dimension dies away, 7^ — > 0, and one recovers the expected bremsstrahlung 
spectrum Uq oc l/fc^. 
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target wavefunction: the latter would be created, or annihilated, by the respective fundamental 
fields only to lowest order in perturbation theory; but at strong coupling this relation could be 
strongly renormalized, in a way which is not under control. It is not clear, though, whether 
such a renormalization will completely wash out the parton picture suggested by perturbation 
theory. For instance, it is known that the entropy of strongly coupled M = A SYM theory at 
finite temperature, which is a direct measure of the number of states, differs only by a factor 
3/4 from the entropy of the corresponding ideal gas, where all the states are associated with 
free fundamental fields. This is suggestive that, at least, the counting of sates is quite similar 
at strong coupling and at weak coupling. 

A similar suggestion emerges from our present results for DIS at strong coupling. When 
analyzed in view of Eq. (U) (with ^ N"^ - I ^ N^, to account for the proper number of 
color degrees of freedom of the fields in = 4 SYM), the estimate ( 4.3^ ) for F2 suggests that, 
for transverse momenta k± ^ Qs(t), the occupation numbers for the fermionic and scalar fields 
saturate at a value of order 1. Note that, for fermions, this value n = 1 is in agreement with the 
maximal occupancy allowed by the Pauli principle. Our analysis gives no direct access to the 
gluon occupancy, but in view of the supersymmetry of the underlying theory, we expect a similar 
saturation value, ~ 1 for k± ^ QsCt), for the gluons as well. This value is quite reasonable: 
at small coupling, one needs a high occupancy ~ l/a^ 3> 1 in order to compensate for 
the weakness of the interactions, but when the coupling is strong, the gluons mutual repulsion 
responsible for saturation becomes important already when ~ 1. 

By pushing this interpretation towards larger virtualities Q"^ S> Q'1{t), where Eq. ( [4.39| ) 
applies, one finds that, when increasing k± above Qs(t), the (scalar and fermion) occupation 
numbers decrease very fast, as n ~ (Q^(r)/A;^)^. Assuming that this behaviour extends to 
gluons as well, we deduce that, at strong coupling, there are essentially no partons at transverse 
momenta above the saturation scale: at large k±, the spectrum converges so fast that any 
interesting convolution involving the integral of n{T,k±) over k± will be dominated by k± ^ 
Qs{t). Equivalently, when the hadron is probed with a given transverse resolution scale Q^, 
almost all partons appear to live at very small values of x, or large values of r = ln(l/2;), above 
the saturation line: r > Ts{Q'^). Since, on the other hand, such small-x partons carry only little 
longitudinal momentum, it is natural to find that the energy-momentum sum rule is dominated 
by those partons living along the saturation line, as shown by Eq. ( [4.37 ). 

It is interesting to notice that the 'phase diagram' in Fig. ^, or, more precisely, a slice of it 
at fixed r, bears some resemblance to a previous scenario for high-energy string scattering in 
fiat space by Amati et al , after replacing the impact parameter b of Ref. by the photon 
virtuality . Specifically, both pictures have in common the dominance of the diffractive over 
the inelastic cross-section at large values of p (with p = In in the present context and p = b 
in the case of Ref. |]l2|) and for sufficiently large values of r. But a closer inspection reveals 
also some important differences between these two pictures, which can be attributed to our use 
of a nontrivial background metric. In particular, and unlike in Ref. in our present analysis 
there is no hierarchy between the elastic and diffractive processes: they all contribute on the 
same footing and simultaneously become large on the 'elastic line' in Eq. ( 4.29| ). 

Consider finally the corner of the 'phase diagram' at ultrahigh energies s > A'^^ which in 
Fig. H is referred to as 'strong gravity'. This is the region where the square of the momentum 
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transferred in the u direction, t = kuk"^ = G^^kuku, which grows with s like tu{s) oc (see 
the Appendix) becomes of the same order as the total energy squared s, meaning that the 
eikonal approximation cannot be trusted anymore. Physically, this is generally interpreted (so 



far, mostly in the context of flat space analyses | p!2| , 47]) as the onset of non-linear gravity 
effects, possibly associated with the formation of a black hole. In this regime, the graviton 
self-interactions cannot be neglected anymore, so in the context of the scattering problem one 
has to resum an infinite number of interacting gravitons in the i-channel. This is reminiscent of 
the non-linear, strong gluon field problem of pQCD in the saturation regime |2^, 0, |8|, 
It would be therefore interesting to see if the similarity between the two phenomena could be 
made even sharper. 
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A. On the breakdown of the eikonal approximation at ultrahigh energies 

In this Appendix, we discuss the breakdown of the eikonal approximation due to large momentum 
transfer at ultrahigh energies. The momentum ku transferred in the u direction can be estimated 
as the M-derivative of the real part of the phase-shift 6 in the S-matrix {S = e**^): k^ ~ 

Consider first the regime at relatively large Q^, such that p = ln(Q^/A^) > 4r/\/A. Then, 
as explained in Sects. 4.2 and 4.3, the phase-shift is controlled by the elementary graviton 
exchange, which yields 

ku-dJ^i)^^4, (A.1) 



(A.2) 



and hence the invariant squared momentum transfer reads 







where we have used G"" = / E? . When this becomes on the order of the total squared energy 
s = {Ff' /v?)s (the rescaling factor E? /v? is necessary to relate the four dimensional energy to 
the ten dimensional one [^]) 

t ^s, (A.3) 

the eikonal approximation breaks down. In the context of DIS, this condition must be considered 
for u ^ Uc = QR^, cf. Eq. ( |3.6D , which is the value that controls the respective cross-section. 
This yields the following critical line 

T = 2p + lnN^, (A.4) 
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above which one expects strong gravity effects, possibly leading to the formation of a black hole. 

Consider similarly the low-Q^ regime, at p < 4r/\/A- Then, the phase shift is associated 
with the single Pomeron exchange, and can be estimated as in Eq. ( 4.26| ) (for both the real and 
the imaginary part): 

5 ~ -^e'^0^^31^ (A.5) 

The u-derivative introduces a factor of order 1/u (since dup = 2/ti), so up to subleading pref- 
actors, one finds 

t ^ G^-{du5f ~ e^-o-^ , (A.6) 

to be compared with s = {E? /u1)s = /E?. This comparison determines a curve 



l- — ]T--^-^ = \nN^, (A.7) 



which smoothly matches the previous line r = 2p + In A'^^, cf. Eq. ( [A.4| ), at a point 

, , / InA-.^ 4inAr4\ 



Together, the two lines in Eqs. ( [A.4| ) and ( [A.§| ) determines the boundary of the applicability of 
the eikonal approximation at ultrahigh energies s > N^, as illustrated in Fig. |3[ 
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